U.S.  DEPARTMENT  OF  COMMERCE 
National  Technical  Information  Service 


AD-A033  821 


THREE-DIMENSIONAL  FRACTURE  ANALYSIS:  PROCEEDINGS 
OF  A WORKSHOP  HELD  AT  BATTELLE'S  COLUMBUS 
LABORATORIES  ON  APRIL  26-28,  1976 


Battelle  Columbus  Laboratories,  Ohio 


30  November  1976 


4DA033821 


r. 


4 


«5  ' 


i > 


; 

t; 

>. 

t 

V 


i 


** 


< 

\ 


t 

i 


i' 

. Unclassified 

1 

, ro 

SECURITY  CLASSIFICATION  OF  THIc  ° AGE  (Wb*n  Data  Knitted) 

REPORT  DOCUMENTATION  PAGE 

READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 

1 

jl.  REPORT  SUMDEH 

ft  F05'R  - TR-  76  - /306. 

2.  GOVT  ACCESSION  NO. 

3.  RECIP'FNT’S  CATALOG  NUMBER 

1 

1 

|4.  1 ITLE  (and  Subtitle) 

5.  TYFE  OF  REPORT  4 PERIOD  COVERED 

i 

Thtee-Dimensional  Fracture  Analysis 

?inal  Report 

April  1,  1976  - Nov.  30,  1976 

« 

6.  PERFORMING  ORG.  REPORT  NUMBER 

f 

7.  AUTHOR/*; 

8.  CONTRACT  OR  GRANT  NUMBER/*; 

L.  E.  Hulbert  (Editor) 

AFOSR-76-294* 

1 

t 

! 

9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

BATTELLE-Columbus  Laboratories 
505  King  Avenue 
Columbus,  Ohio  4320'J 

10.  PPOGPAM  F.LEMENT.  PROJECT,  TASK 
APEA  A WORK  UNIT  NUMBERS 

681307  2307B1 

61102F 

• 

11.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH/NA 
BLDG  410 

BOLLING  AIR  FORCE  BASE,  D C 20332 

12.  REPORT  DATE 

November  30,  1976 

> 

13.  NUMBER  OF  PAGES 
116 

• 

14  MONITORING  AGENCY  NAME  & AODRESSf/4  dillerent  Irom  Controlling  Oil  ice) 

IS.  SECURITY  CLASS,  (ol  th‘e  report) 

Unclassified 

15a.  DECLASSIFICATION/ DOWN  GRADING 
SCHEDULE 

i 

\ 

16.  DISTRIBUTION  STATEMENT  (ol  this  Report) 

i 

Approved  for  public  release.  Distribution  unlimited. 

V 

1 

< 

S 

j 

17.  DISTRIBUTION  STATEMENT  (of  the  abstract  entered  in  Block  20,  if  different  fret 

n Report) 

i 

< 

i 

j 

1 

18.  SUPPLEMENTARY  NOTES 

,1 

1 

T 

* 

f 

I 19.  KEY  WORDS  (Continue  on  reverse  side  if  necessary  and  identify  by  block  number) 

5 

racture  Workshop  Quarter  point  Cracked  Elements 

3-d  Fracture  Crack  Tip  Phenomena 

Fracture  Analysis  Government  Fracture  Problems 

y rid  (.racked  Elements  Benchmark  Fracture  Problems 

Boundary  Integral  Analysis 

i 

I 20.  ABSTRACT  (Continue  on  reverse  aide  If  necessary  ana  identify  by  block  number^ 

A three  day  workshop  on  analysis  of  three-dimensional  stress  states  around 
cracks  was  held  at  Battelle's  Columbus  Laboratories  on  April  24-26  1976. 

This  workshop  was  attended  by  outstanding  American  scientists  and  government 
representatives  concerned  with  three-dimensional  analysis  of  fracture.  An 
intense e effort  was  made  to  identify  and  invite  all  of  the  active  researchers 
m this  field.  The  result  of  having  virtually  all  of  the  known  researchers 

— Vt»  - ..., 


J 


Unclassified 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fKh7n  Irate  tntrredi 


DD  | JAN^S  1473  EDITION  OF  I NOV  65  IS  OBSOLETE 


Unclassified — i— 

SECORITl  CLASSIFICATION  OF  THIS  PACECWien  Date  Entarad) 


attend  the  workshop  presented  a unique  opportunity  to  establish  the 
state-of-the-art  in  this  field. 

The  current  report  includes  abstracts  of  the  papers  presented,  short 
summaries  of  the  discussions,  summaries  of  remarks  by  government 
representatives  concerning  their  fracture  problems,  and  a definition 
of  three  benchmark  problems  which  are  to  be  used  as  standards  of 
comparisons  for  future  analysis  methods. 

This  workshop  was  jointly  supported  by  the  Air  Force  Office  of  Scientific 
Research  and  the  Transportation  Systems  Center,  Department  of  Transportatior , 
with  contributing  support  of  the  Army  Research  Office. 


Unclassif ied 

SECURITY  CLASSIFICATION  OF  Tf.lS  PAG  E(W*on  t nt*red) 


..  v's 


Ill 


£ AC’tiSIMJtr 

. Wktta  Sidln 
i lafl  SKtlM  □ 

I '•'.r'-fMCEl  □ 

j ;■*:?:«!»» 

i , 

'<Kt<  

• JTfiiwriOK/msLAiiun  wies 

}'  i>Yf  AVAIL  ind/ar  VECiAi" 


\ 


THREE-DIMENSIONAL  FRACTURE  ANALYSIS 


Edited  By 

L.  E.  Hulbert 
Battelle-Columbus 


November  30,  1976 


PROCEEDINGS  OF  A WORKSHOP 
HELD  AT 

BATTELLE'S  COLUMBUS  LABORATORIES 
April  26-28,  1976 


For 

DEPARTMENT  OF  TRANSPORTATION  (TSC) 
U.  S,  AIR  FORCE 

U.  S.  ARMY  RESEARCH  OFFICE,  DURHAM 


D D C 

H35M3E  nrpi 


fEISMITElUi 

D 


DISTRIBUTION  STATEMENT  _A_ 
Approved  loi  public*  leloase; 
Distiibution  Unliiiiited 


«*ra*r 


•Jet* 


I 

I 


LIST  OF  PARTICIPANTS  

INTRODUCTION  

SUMMARY  AND  CONCLUSIONS 

ACKNOWLEDGEMENTS  

SESSION  I:  FINITE  ELEMENT  METHODS 

Crack  Tip  Plasticity  Calculations  With  Quadratic 

"Quarter  Point"  Finite  Elements,  By  S.  E.  Benzley  . , 

Three-Dimensional  Finite  Element  Analysis  for  the 

Through  Crack  Problem,  By  P.  D.  Hilton  

Quarter  Point  Singularity  Elements  and  Their 
Application  to  Three-Dimensional  Fracture  Analysis, 

By  R.  S.  Barsoum  . , 

On  a 3-D  "Singularity-Element"  for  Computation  of 

Mixed-Mode  Stress  Intensities,  By  3.  N.  Atluri,  et  al  

Comments  on  3D  Finite  Element  Elastic  Crack  Analysis, 

By  D.  M.  Tracey  

An  Evaluation  of  the  Quadratic  Isoparametric  Singularity 
Element,  By  I.  M.  Bloom 

On  Some  Fracture  Analysis  Results  at  the  Berkeley 

Nuclear  Laboratories,  By  C.  H.  A.  Tow-.’ey 

SESSION  II:  CRACK  TIP  PHENOMENA 

Stress  Intensity  Estimates  for  Three-Dimensional  Cracked 
Body  Problems  by  the  Frozen  Stress  Photoelastic  Method, 

By  C.  W.  Smith  

Observations  of  Crack  Tip  Processes,  By  G.  T.  Hahn  

Use  of  Cyclic  Growth  Tests  to  Infer  Stress  Intensity, 

By  J.  E.  Collipriesi  

Comments  for  Workshop  on  Three-Dimensional  Fracture 

Analysis,  J.  L.  Swadlow . 

Crack  Tip  Fields  in  Steady  Crack  Growth  with  Strain 

Hardening,  By  J.  W.  Hutchinson - . . . 


Page 

4 

1 

1 

2 

3 

6 

9 

11 

14 

17 

22 

d5 

34 

42 

45 

48 


TABLE  OF  CONTENTS 
(Continued) 


\/ 


Page 

Some  Properties  of  Finite  Element  Approxi?.  ations  of 
Elliptic  Problems  on  Domains  with  Cracks  and  Corners, 

By  J.  T.  Oden 49 

Near  Field  Behavior  and  Crack  Growth,  By  G.  C.  Sih 50 

Stress  Intensity  Factor  Measurements  for  Corner  Cracked 

Holes,  By  A.  F.  Grandt,  Jr 51 

SESSION  III:  GLOBAL  FUNCTION  METHODS 

Development  of  Procedures  for  Analyzing  Stresses  in 

Cracked  Bodies  of  Various  Shapes,  By  J.  C.  Bell  . 52 

Boundary- Integral  Equation  Analysis  of  Surface 

Cracks,  By  T.  A.  Cruse 57 

Subsurface  Elliptical  Flaws,  By.  A.  S.  Kobayashi,  et  al  65 

Stress  Intensity  Factors  for  a Pressurized  Thick-Wall 
Cylinder  with  a Part-Through  Circular  Surface  Flaw  - 

Compliance  Calibration  and  Collocation  Method,  By  M.  A.  Hussain  . . 68 

On  The  Three-Dimensional  Theory  of  Fracture,  By  E.  S.  Folias  ...  70 

Some  Unsolved  Singularity  Problems,  By  M.  L.  Williams  72 

SESSION  IV:  FINITE  ELEMENT  AND  FINITE  DIFFERENCE  METHODS 

Fundamental  Study  of  Crack  Initiation  and  Propagation, 

By  M.  L.  Wilkins 73 

Application  of  an  Influence  Function  Method  for  Three- 

Dimensional  Elastic  Analysis  of  Cracks,  By  P.  M.  Besuner 74 

Elliptical  Crack,  Surface  Flaw,  and  Flaws  at  Fastener 

Holes,  By  R.  C.  Shah 77 

An  Assessment  of  Near  Tip  Modeling  for  2D  and  3D  Crack 

Problems,  By  C.  F.  Shih 81 

Hybrid  Models  for  Three-Dimensional  Crack  Elements, 

By  T.  H.  H.  Pian 82 

The  Finite  Element  Alternating  Method  for  Analysis  of 

Complex  Three-Dimensional  Crack  Problems,  By  F.  W.  Smith  84 


V* 


TABLE  OF  CONTENTS 
(Continued) 


SESSION  V:  PRACTICAL  PROBLEMS  OF  INTEREST  TO  GOVERNMENT  AGENCIES 


U.  S.  Army  (Watervliet  Arsenal) , By  M.  A.  Hussain 


U.  S.  Army  (AMMRC) , By  D.  M.  Tracey 


U.  S.  Air  Force  (AFFDL) , By  H.  A.  Wood 
U.  S.  Air  Force  (AFML) , By  T-  Nicholas  . 
U.  S.  Air  Force  (AFRPL) , By  R.  Peeters  . 
U.  S.  Air  Force  (AFOSR) , By  W.  J.  Walker 


NASA  (Lewis  Research  Center),  By  C.  C.  Chamis 


Department  of  Transportation  (TSC) , By  D.  McConnell,  E.  Savage  . . 91 


Oak  Ridge  National  Laboratory  (ERDA),  By  G.  Smith 
SESSION  VI:  BENCHMARK  PROBLEMS 


SUMMARY 100 


CRITERIA ^ 101 

BENCHMARK  NUMBER  1 THE  SURFACE  FLAW 102 

BENCHMARK  NUMBER  2 THE  CORNER  CRACKED  HOLE  105 

BENCHMARK  NUMBER  3 THE  COMPACT  SPECIMEN  107 


\-f . 'v' 


LIST  OF  PARTICIPANTS 


S.  N.  Atluri 

R.  S.  Barsoura 
J.  C.  Bell 

S.  E.  Benzley 
P.  M.  Besuner 
J.  M.  Bloom 

D.  Broek 

C.  C.  Chamis 
J.  E.  Collipriest 

T.  A.  Cruse 
R.  M.  Engle 

E.  S.  Folias 
A.  F.  Grandt 
G.  T.  Hahn 

C.  Hayes 

P.  D.  Hilton 
A.  T.  Hopper 

L.  E.  Hulbert 

J.  W.  Hutchinson 

M.  A.  Hussain 
T.  J.  Johns 
R.  E.  Jones 

M.  F.  Kanninen 
A.  S.  Kobayashi 
G.  E.  Maddux 

D.  McConnell 
T.  Nicholas 
J . T . Oden 


Organization 

Georgia  Institute  of  Technology 
Combustion  Engineering,  Inc. 
Battelle-Columbus  Laboratories 
Sandia  Laboratories 
Fracture  Analysis  Associates 
Babcock  and  Wilcox 
Battelle-Columbus  Laboratories 
NASA,  Lewis  Research  Center 
Rockwell  International 
Pratt  and  Whitney  Aircraft,  Inc. 

U.  S.  A.  F.  Flight  Dynamics  Laboratory 

The  University  of  Utah 

U.  S.  A.  F.  Materials  Laboratory 

Battelle-Columbus  Laboratories 

U.  S.  A.  F.  Office  of  Scientific  Research 

Lehigh  University 

Battelle-Columbus  Laboratories 

Battelle-Columbus  Laboratories 

Harvard  University 

U.  S.  Army  Watervliet  Arsenal 

Battelle-Columbus  Laboratories 

Boeing  Aerospace  Corporation 

Battelle-Columbus  Laboratories 

University  of  Washington 

U.  S.  A.  F.  Flight  Dynamics  Laboratory 

Transportation  Systems  Center,  DOT 

U.  S.  A.  F.  Materials  Laboratory 

University  of  Texas  at  Austin 


Mil 


\atxtti-xsa -evSiZ 


LIST  OF  PARTICIPANTS 
(Continued) 


Organization 


R. 

Peeters 

U.  S.  A.  F.  Rocket  Propulsion  Laboratory 

T. 

H.  H.  Pian 

Massachusetts  of  Technology 

; 

C. 

H,  Popelar 

The  Ohio  State  University 

: 

I. 

S.  Raju 

NASA,  Langley  Research  Center 

E. 

F.  Rybicki 

Battelle-Columbus  Laboratories 

;■ 

S. 

G.  Sampath 

Battelle-Columbus  Laboratories 

‘ 

E. 

Savage 

Transportation  Systems  Center,  DOT 

R. 

C.  Shah 

Boeing  Aerospace  Company 

j 

C. 

F.  Shin 

General  Electric  Company 

s 

G. 

C.  Sih 

Lehigh  University 

< 

C. 

W.  Smith 

Virginia  Polytechnic  Institute  and 

J 

State  University 

] 

j 

F. 

W.  Smith 

Colorado  State  University 

1 

J 

G. 

Smith 

Oak  Ridge  National  Laboratory 

i 

/ 

i 

J. 

L.  Swedlow 

Carnegie-Mellon  University 

i 

P. 

Tong 

Transportation  Systems  Center,  DOT 

C. 

H,  A.  Townley 

C.  E.  G.  B.  Berkeley  Nuclear  Laboratories 

* 

(England) 

D. 

M.  Tracey 

U.  S.  Army  Materials  and  Mechanics 

* 

Research  Center 

| 

M. 

L.  Williams 

University  of  Pittsburgh 

X 

W. 

J.  Walker 

U.  S.  A.  F.  Office  of  Scientific 

’ k- 

Research 

\ 

M. 

L.  Wilkins 

Lawrence  Livermore  Laboratory 

s 

H. 

A . Wood 

U.  S.  A.  F.  Flight  Dynamics  Laboratorv 

* 

THREE-DIMENSIONAL  FRACTURE  ANALYSIS 


Edited  by 


L.  E.  Hulbert 


INTRODUCTION 


A three  day  workshop  was  held  at  Battelle's  Columbus  Laboratories 
on  April  26-28,  with  the  purpose  of  bringing  together  the  nation's  leading 
investigators  of  three-dimensional  fracture  to  discuss  togetb  r the  pre- 
sent status  of  three-dimentional  fracture  analysis.  The  Workshop  was  orga- 
nized and  conducted  by  Battelle  under  the  joint  sponsorship  of  the  U.  S. 

Air  Force  and  the  Department  of  Transportation,  with  contributing  support 
of  the  U.  S.  Army  Research  Office. 

The  Workshop  was  organized  in  six  sessions  with  the  first  four 
devoted  to  technical  papers  and  discussions,  one  session  in  which  govern- 
ment representatives  discussed  their  problems  and  one  session  devoted  to 
the  definition  of  benchmark  problems. 

This  report  presents  a summary  of  the  Workshop.  The  four  paper 
sessions  are  presented  in  the  form  of  extended  abstracts  written  by  each 
author  with  summaries  of  discussions  after  each  abstract.  The  remarks  of 
the  government  representatives  are  summarized  together  with  discussions  of 
each  presentation.  A summary  is  also  given  of  the  session  in  which  stan- 
dards for  benchmark  problems  were  discussed  and  three  benchmark  problems 
were  chosen.  These  problems  are  described  in  detail. 


SUMMARY  AND  CONCLUSIONS 

In  organizing  chis  Workshop,  Battelle  made  an  intensive  effort 
to  identify  all  of  the  scientists  (exclusive  of  graduate  students)  in  the 
United  States  who  had  been  conducting  significant  research  into  the  analysis 
of  thrc e-dimensional  stress  states  around  cracks.  Although  it  is  possible 
that  some  such  scientists  were  not  identified,  no  one  nas  been  brought  to 
our  attention  as  of  the  date  of  th>'s  report.  Remarkably,  all  scientists 
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identified  as  working  on  three-dimensional  aspects  of  fracture  accepted 
Battelle’s  invitation  and  attended  the  Workshop. 

A similar  effort  was  made  to  identify  and  invite  those  govern- 
ment representatives  interested  in  the  subject.  Because  of  unavoidable  con- 
flicts, four  of  these  representatives  could  not  attend.  However,  nearly  all 
of  the  govenment  agencies  interested  in  fracture  were  represented.  Thus, 
this  Workshop  represented  a unique  event  in  the  history  of  fracture  analysis 
in  terms  of  defining  precisely  the  current  state  of  three-dimensional  frac- 
ture analysis  and  defining  future  directions  for  research. 
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FINITE  ELEMFNT  ANALYSIS  METHODS 
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CRACK  TIP  PLASTICITY  CALCULATIONS  WITH  QUADRATIC 
"QUARTER  POINT"  FINITE  ELEMENTS 

by 

S.  E.  Benzley 

Sandia  Laboratories 
Albuquerque,  New  Mexico 


The  dominant  plastic  singularity  at  the  tip  of  a crack  for  a pure 
power  law  hardening  material  has  been  expressed  as: 
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n + 1 
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where  (r,  0)  are  polar  coordinates  centered  at  the  crack  tip,  and  e are 
the  near  field  stress  and  strain  respectively,  n is  the  hardening  exponent 
chosen  to  represent  the  experimentally  determined  stress-stiain  curve,  and 

A A 

a, , and  e are  functions  giving  the  0 dependence.  The  explicit  determination 

A A 

of  the  a^(0)  and  (0)  terms  involves  the  solution  of  ? nonlinear,  fourth 
order  ordinary  differential  equation  for  single  mode  problem;  in  plain  stress 
or  plane  strain.  Mixed  mode  solutions  become  even  more  complex. 

Efforts  to  characterize  this  singularity  with  finite  elements  have 
involved  (1)  defining  the  region  around  the  crack  with  a high  concentration  of 
finite  elements  or  (2)  developing  special  crack  tip  elements  that  incorporate 
the  singularity  of  Equation  (1).  Boi_h  "displacement"  and  "hybrid"  formulations 
have  been  successful. 

The  quadratic  isoparametric  finite  element  has  proven  itself  very 
valuable  in  two  dimensional  fracture  mechanics  problems.  This  is  because 
a — singularity  can  be  established  in  the  strain  displacement  matrix  by 


moving  selected  midside  nodes  to  the  element  "quarter  points".  This  idea 
has  been  extended  to  three-dimensional  problems  using  quadratic  brick  and 
wedge  type  elements  and  again  moving  midside  nodes  to  quadratic  point 
positions.  Reported  here  is  the  use  of  this  technique  on  elastic-plastic 
problems. 


Considering  the  "strain  energy"  (i.e.  / dv)  the  crack 

tip  and  the  near  field  representation  of  stress  and  strain  as  given  by 
Equation  1,  one  may  conclude  that  the  strain  energy  for  plastic  (as  well  as 
elastic)  behavior  has  1/r  characteristic  singularity  at  a crack  tip,  i.e., 


S.E. 


(0)dv 


(2) 


The  stiffness  matrix,  K,  used  in  finite  element  calculations  is 
formed  from  the  strain  energy  terms  and  may  be  expressed  as 


K = / BTDB  dv  , (3) 


where  B is  the  matrix  transforming  displacements  to  strains  and  I)  is  a 
material  matrix.  The  quarter  point  adjustment  has  the  effect  of  establishing 
a 1/r  singularity  in  the  B matrix,  thus 

K = / - BTDB  dv  . (4) 

r 

The  above  heuristic  argument  thus  establishes  the  quarter  point  element  as  a 
possible  method  for  nonlinear  crack  tip  computations  in  both  two  and  three 
dimensions. 

Comparisons  have  been  made  on  a two  dimensional  problem  using  this 
technique  and  a solution  incorporating  a special  plastic  singular  element. 
This  comparison  indicates  that,  for  the  case  of  two  dimensional  single  mode 
problems,  the  quarter  point  calculations  give  very  reasonable  results. 

As  previously  stated,  this  approach  to  crack  tip  plasticity  is 
directly  applicable  to  three-dimensional  calculations.  In  fact,  any  three- 
dimensional  elasto-plastic  computer  program  that  incorporates  a quadratic 
element  can  be  used  as  it  stands,  to  solve  crack  tip  plasticity  problems. 
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Discussion 

P.  D.  Hilton  pointed  out  that  agreement  between  the  results  pre- 
sented and  those  of  reference  solutions  was  good  on  a gross  behavior  scale, 
but  asked  whether  any  comparisons  on  a more  local  scale,  such  as  a J-integral, 
had  been  made.  S.  E.  Ber.zley  replied  that  to  date  those  types  of  comparisons 
had  not  been  made. 

C.  C.  Chamis  questioned  the  meaning  of  "engineering  solution"  as 
described  in  the  presentation.  Chamis  contended  that  engineering  solutions 
should  be  compared  with  experimental  data  rather  than  other  numerical  solu- 
tions. 

R.  S.  Barsoun  asked  what  results  had  been  obtained  for  cases  with 
flaw  curvature.  Benzley  replied  that  limited  work  on  that  problem  did  not 
produce  the  good  correlation  that  was  presented  and  results  near  the  free  edge 
gave  the  least  correlation. 

J.  Swedlow  pointed  out  that  crude  analyses  are  sensitive  to  any 
curvature  of  the  three-dimensional  crack  surface. 

S.  G.  Sampath  inquired  if  interelement  compatibility  was  satisfied 
near  the  crack  tip.  Benrley  replied  that  compatible  elements  were  used  in 
all  of  the  results  presented.  Chamis  inquired  how  equilibrium  was  satisfied. 
Benzley  replied  equilibrium  was  satisfied  approximately  by  the  finite  element 
method. 


THREE-DIMENSIONAL  FINITE  ELEMENT  ANALYSIS 
FOR  THE  THROUGH  CRACK  PROBLEM 


P.  D.  Hilton 

Institute  of  Fracture  and  Solid  Mechanics 
Lehigh  University 
Bethlehem,  Pennsylvania 


The  problem  of  a through  crack  in  a plate  subjected  to  in~plane 
loading  has  traditionally  been  treated  using  the  two-dimensional  idealiza- 
tions of  plane  stress  or  plane  strain.  When  significant  yielding  occurs, 
the  observed  response  exhibits  some  three-dimensional  characteristics  in 
conjunction  with  stable  crack  growth,  i.e.,  crack  front  curvature  and  the 
development  of  shear  lips.  Our  long  range  goal  is  to  model  this  problem 
and  analyze  these  effects.  Some  initial  work  in  the  area  of  three-dimen- 
sional elastic  finite  element  analysis  has  been  performed.  The  three- 
dimensional  elastic-plastic  effort  is  underway. 

Elastic  Singularity  Elements 


Based  on  previous  finite  element  accomplishments  in  two-dimen- 
sional fracture  mechanics,  a number  of  approaches  for  the  development  of 
three-dimensional  singular  elements  for  crack  problems  are  available.  We 
have  chosen  to  incorporate  the  concepts  utilized  in  Benzley's  "enriched" 
elements*  as  the  basis  of  a three-dimensional  singular  element.  These 
elements  are  characterized  by  an  assumed  form  for  the  dispx.. cement  field 
which  contains  the  first  term  of  the  asymptotic  solution  (corresponding  to 
the  square  root  singularity  in  the  near  tip  strain  field)  and  the  stancard 
polynomial  terms  associateu  with  Isoparametric  elements.  The  stress  intensity 
factor (s)  is  treated  as  an  additional  unknown  (generalized  nodal  displacement 

* S.  E.  Benzley,  "Representation  of  Singularities  with  Isoparametric 
Finite  Elements",  Int.  J.  Num.  Meth.  Eng. , Vol.  8,  pp  537—54 5 (1974). 
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component)  to  be  determined  by  the  finite  element  analysis.  For  the  three- 
dimensional  analog,  variation  of  the  stress  intensity  factor (3)  along  the 
crack  front  is  approximated  in  a piece-wise  fashion  by  polynomials,  i.e.,  a 
polynomial  distribution  is  assumed  over  the  edge  of  each  singular  element 
which  coincides  with  a portion  of  the  crack  front.  The  values  of  the  stress 
intensity  factor (s)  at  nodes  along  the  crack  fron  are  unknowns  determined  by 
the  minimization  of  the  potential  energy  (as  approximated  by  the  finite 
element  method) . 

The  use  of  these  enriched  elements  which  contain  both  the  singu- 
larity solution  and  the  nonsingular  polynomial  approximations  for  the  dis- 
placement field  is  particularly  suited  for  application  to  three-dimensional 
crack  problems.  These  elements  are  able  to  model  both  the  interior  region, 
where  the  singular  solution  dominates  in  the  vicinity  of  the  crack  edge,  and 
the  region  near  the  free-surface  in  which  the  stress  intensity  factor  may  be 
small  or  zero  and  the  large  (but  nonsingular)  out-of-plane  shearing  stresses 
become  important. 

Initial  results  indicate  that  the  predictions  for  the  variation  of 
the  stress  intensity  factor  along  the  crack  front  can  be  sensitive  to  the 
radius  of  the  singular  elements  and  that  corresponding  care  must  be  taken  i.i 
the  development  of  grid  patterns  to  obtain  consistent  values.  Calculations 
have  been  performed  which  demonstrate  the  influence  of  plate  thickness  and  of 
crack  front  geometry  on  the  stress  intensity  factor  distributions. 

Preliminary  Elastic-Plastic  Calculations 

A three-dimensional,  incremental  theory  plasticity,  finite  element 
computer  program  has  been  developed  to  study  influences  of  plasticity  for  the 
through  crack  problem  (At  present,  no  special  treatment  of  the  crack  edge 
behavior  is  incorporated.).  Preliminary  results  indicate  that  the  tensile 
stress  (0^)  ahead  of  the  crack  edge  varies  more  rapidly  along  that  edge  when 
plasticity  is  accounted  for.  In  particular,  elastic-plastic  predictions  for 
0^  are  larger  than  elastic  predictions  in  the  vicinity  of  the  mid-plane  of  the 
plate  and  smaller  than  elastic  predictions  near  the  free  surface.  This  result 
is  consistent  with  experimental  observations  on  crack  front  curvature. 


Further  calculations  modeling  increments  of  crack  front  growth  are 
planned  in  an  effort  to  understand  the  conditions  for  stability  of  this 
growth  and  the  corresponding  development  of  both  crack  front  curvature  and 
shear  lips. 

Discussion 


Question:  Why  was  at  the  free  surface  set  to  zero  in  the 

finite  element  'alculations?  Motivation  for  doing  this  was  questioned  on 
the  basis  that  K^.  at  the  free  surface  is  not  known.  Hilton  showed  results 
indicating  that  when  at  the  free  surface  was  not  present,  the  resulting 
Kj  at  the  center  changed  very  little.  Another  discussion  topic  was  that  the 
three-dimensional  result  was  10  percent  higher  at  the  center  of  the  plate 
than  the  result  for  the  two-dimensional  case.  The  effects  of  crack  front 
curvature  and  tunneling  were  brought  up  as  being  important  in  affecting  the 
distribution  of  values.  A question  was  asked  concerning  attempts  to  find 
a crack  front  shape  that  produced  a constant  value  of  K through  the  thickness. 
Hilton  answered  that  preliminary  attempts  to  do  this  have  not  been  successful. 
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QUARTER  POINT  SINGULARITY  ELEMENTS  AND  THEIR 
APPLICATION  TO  THREE-DTMENSIONaL  FRACTURE  ANALYSIS 


Roshdy  S.  Barsoum 

Combustion  Engineering,  Inc. 
Windsor,  Connecticut 


Special  quarter  point  singularity  finite  elements  are  applied  to 
the  elastic  fracture  analysis  of  three-dimensional  surface  and  through  cracks. 
The  applicability  of  the  elements  to  the  elastic-plastic  fracture  analysis 
is  also  discussed, 

(1  2) 

Recently,  it  was  discovered  by  the  author  * and  independently  by 
(3) 

Kenshell  and  Shaw  that  second  order  isoparametric  element  formulation 
allows  a 1/^r  strain  singularity  when  the  side  nodes  are  placed  at  the  1/4 
point  position.  This  singularity  is  exactly  that  of  the  elastic  crack  solu- 
tion. Due  to  the  existence  of  constant  strain  terms,  rigid  body  motion  and 
compatibility  with  surrounding  elements,  they  converge  to  the  exact  solution 
and  admit  easily  loadings  such  as  thermal  loadings.  It  was  shown  that  ex- 
tremely high  accuracy  in  stress  prediction  close  to  the  crack  tip  can  be 
obtained  using  a reasonable  size  mesh  (4).  The  same  theory  was  also  applied 
to  the  8-noded  isoparametric  thick  shell  element  in  the  case  of  a through 
crack  in  a plate  subject  to  bending^.  The  comparison  of  this  analysis  with 
three-dimensional  analytical  solutions ^ lead  to  extremely  high  accuracy 
since  the  formulation  of  these  elements  is  based  on  Reissner's  plate  theory. 

Quarter  point  three-dimensional  elements  were  applied  to  the 
analysis  of  a semi-elliptical  crack  subject  to  thermal  shock. ^ There  is  no 
analytical  solution  for  the  thermal  shock  problem,  however,  comparisons  with 
results  obtained  by  the  alternating  method  technique  for  the  case  of  axial 
loads  were  in  good  agreement. 

In  evaluating  the  accuracy  of  the  results  of  the  quarter  point 
elements,  it  was  found  that  the  triangular  form  cf  these  elements  leads  to 
-• — 1 *• - -1 — *-u"  1 --  >8)  'fhis  led  to  further 


far  superior  results  than  the  rectangular  form 


(8  9)  (9) 

investigation  of  the  element  singularity  ’ , It  was  found  that  the 

triangular  form  of  the  elements  contain  both  l//r  singularity  and  1/r 
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singularity.  The  1/Jr  singularity  is  achieved  when  the  crack  tip  is  not 
allowed  to  blunt,  while  the  1/r  singularity  is  achieved  when  the  crack  tip 

is  allowed  to  blunt  by  having  multiple  independent  nodes  at  the  tip.  Tests 
on  power  law  Hardening  materials  showed  that  the  stress  distribution  at  the 
integration  points  is  in  good  agreement  with  other  results  obtained  using 
special  power  law  crack  tip  elements The  only  problem  remaining  with 
the  three—. imensional  elastic-plastic  analysis  of  cracks  is  the  large  computer 
time  required.  This  is  due  to  the  fact  that  fracture  analysis  requires  a 
rather  refined  mesh  with  large  number  of  degrees  of  freedom  in  order  to  ob- 
tain a meaningful  answer.  At  the  moment  this  problem  remains  to  be  resolved. 
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ON  A 3-D  "SINGULARITY- ELEMENT"  FOR  COMPUTATION  OF 
MIXED-MODE  STRESS  INTENSITIES 


Satya  N.  Atluri,  K.  Kathiresan,  and  M.  Nakagaki 

Georgia  Institute  of  Technology 
Atlanta,  Georgia 


Two  topics  will  be  briefly  discussed  in  the  presentation.  The 
first  deals  with  a finite-element  procedure  for  the  calculation  of  combined 
modes  I,  II,  and  III  stress  intensity  factors,  which  vary,  along  an  arbi- 
trarily curved  three-dimensional  crack  front  in  a structural  component.  The 
finite-element  model  is  based  on  a modified  variational  principle  of  potential 
energy  with  relaxed  continuity  requirements  for  displacements  at  the  inter- 
element boundary.  The  variational  principle  is  a three-field  principle, 
with  the  arbitrary  interior  displacements  for  the  element,  interelement 
boundary  displacements,  and  element  boundary  tractions  as  variables.  The 
unknowns  in  the  final  algebraic  system  of  equations,  In  the  present  displace- 
ment hybrid  finite-element  model,  are  the  nodal  displacements  and  the  three 
elastic  stress-intensity  factors  at  nodes  along  the  crack  front.  Inter- 
element displacement  compatibility  is  satisfied,  by  assuming  an  indenpendent 
interelement  boundary  displacement  field,  and  using  a Lagrange  Multiplier 
technique  to  enforce  such  interelement  compatibility.  These  Lagrange  Multi- 
pliers, which  are  physically  the  boundary  tractions,  are  assumed  from  an 
equilibrated  stress  field  derived  from  three-dimensional  Beltrami  (or  Maxwell- 
Morera)  stress  functions  that  are  complete.  However,  considerable  care 
should  be  exercised  in  the  use  of  these  stress  functions  such  that  the  stresses 
produced  by  any  of  these  stress  function  components  are  not  linearly  dependent. 
Since  the  method  is  based  on  a rigorous  variational  principle,  which  enforces 
at  least  on  an  average  the  conditions  of  interelement  displacement  and  trac- 
tion continuity  when  Jr  type  displacements  are  included  in  the  near-tip  region, 
the  convergence  of  the  finite-element  solution  for  nodal  displacements  as 
well  as  the  stress- intensity  factors  is  established  mathematically.  The 
geometry  of  the  "basic  element"  used  presently,  is  a 20  node  isoparametric 
"brick"  element,  with  60  degrees  of  freedom  per  element.  The  relevant 
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matrices  are  evaluated  numerically,  using  non-product  type  quadrature 
formulae  with  proper  mathematical  transformations  being  used  when  singular- 
type  functions  are  encountered  in  stresses  and  strains  in  the  near-tip 
region. 

The  ease  with  which  the  above  'singular’  element  can  be  implemented 
in  existing  general  purpose,  efficient,  3-D  finite  element  codes  will  be 
discussed. 

The  second  topic  deals  with  research  in  progress  on  the  applica- 
tion of  the  "edge-function"  method  to  3-D  crack  problems.  The.  edge-function 
method  may  be  described  as  a piecing  togethei  of  "asymptotic"  solutions  to 
the  governing  differential  equations  for  the  several  parts  of  a domain  to 
satisfy  the  boundary  conditions  in  a discrete  least  squares  sense.  The 
edge-function  method  for  the  present  problem  leads  to  •.  'super-element'  with 
embedded  singularities  similar  to  the  one  described  above. 

Some  sample  problems  to  check  the  developed  proced  ires  will  be  dis- 
cussed. One  of  these  is  the  problem  of  a through  the  thickness  crack  in  a 
finite-width  plate  subjected  to  combined  tension,  in-plane  shear  and  out-of- 
plane shear.  In  this  problem,  the  boundary  layer  effects  at  the  intersection 
of  the  through  crack  with  the  surface  of  the  plate  will  be  discussed.  The 
second  problem  is  that  of  a sandwich  plate  with  a part-through  and  a debonding 
crack.  It  is  assumed  that  two  outer  layers  are  bonded  through  an  adhesive 
of  constant  thickness,  and  one  of  the  outerplates  is  assumed  to  contain  a 
through  crack  of  finite  length.  The  results  obtained  for  this  problem  are 
compared  with  the  two-dimensional  analysis  results  of  Erdogan  and  Arin.  An- 
other set  of  problems  to  be  discussed  is  that  of  circular  or  el1  ptical  shaped 
cracks  that  are  either  embedded  or  at  the  corner  of  a tension  bar  of  square 
cross  section.  The  obtained  results  are  compared  with  those  of  Cruse  using 
a boundary- integral  method. 
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DISCUSSION 


The  discussion  following  this  presentation  focused  on  the  question 
of  whether  the  three  stress  intensity  factors  are  sufficient  to  explain  all 
of  the  crack  tip  behavior  and  the  dependence  of  K values  the  technique 
used  to  evaluate  K.  It  was  pointed  out  by  G.  Sih  that  results  at  the  free 
surface  suggest  the  K^,  K, ^ , and  K^^  may  not  be  enough  to  describe  the 
behavior  and  this  question  should  not  be  dismissed. 
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COMMENTS  ON  3D  FINITE- ELEMENT  ELASTIC  CRACK  ANALYSIS 


by 

Dennis  M.  Tracey 

Array  Materials  and  Mechanics  Research  Center 
Watertown,  Massachusetts 


The  discussion  concerns  the  writer's  experience  in  adapting  the 
finite-element  method  to  3D  elastic  crack  analysis.  The  work  described  was 
performed  to  satisfy  the  need  for  an  accurate  means  to  apply  LEFM  concepts 
to  complex  structures,  and  also  to  provide  a way  to  accurately  evaluate 
specimens  used  for  crack  growth  and  fracture  testing  purposes.  The  choice 
of  the  finite-element  method  to  address  these  problems  naturally  stems  from 
the  bacic  versatility  of  the  method:  the  most  complex  loading  involving  the 

most  awkward  geometry  can  be  modeled  in  a systematic,  standardized  fashion. 

We  will  be  concerned  here  with  those  features  of  a finite-element  model  which 
influence  the  accuracy  of  crack  solutions.  Primary  attention  is  gi^en  to 
the  proper  choice  of  deformation  modes  of  elements  used  at  the  crack  front, 
recognizing  the  singular  nature  of  the  deformation  there.  The  necessity  of 
precision  in  modeling  the  actual  hardware  load  state  and  constitutive  behavior 
is  emphasized,  as  this  dictates  the  degree  to  which  the  advances  in  finite- 
e’ement  formulation  are  practically  beneficial. 

Analytical  studies  have  revealed  that  the  inverse  square  root  singu- 
larity governs  at  a 3D  crack  front,  and,  excluding  a Mode  III  crack  face 
sliding  situation,  in  an  asymptotic  sense  a plane  strain  deformation  state 
exists  at  the  front.  (One  important  qualification  to  this  is  at  the  point 
where  a crack  front  intersects  a boundary,  firm  knowledge  of  the  singularity 
there  is  lacking.)  Experience  with  element  assumed  fields  of  the  standard 
polynomial  type  has  supported  what  theoretically  has  been  predicted:  con- 

vergence to  an  accurate  singularity  solution  cannot  be  achieved  with  nen- 
singular  deformation  representations.  This  experience  has  led  to  the  invention 
of  numerous  techniques  for  deduction  of  the  stress  intensity  factor (s)  from 
fundamentally  inaccurate  crack  solutions.  The  discucr^on  here  will  not  elab- 
orate on  these  techniques,  but  instead,  will  address  the  use  of  singularity 
elements  to  achieve  accurate  crack  solutions,  from  which  unambiguously  follows 
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the  stress  intensity  factor  distribution  or  any  other  field  quantity. 

While  a plane  strain  square  root  singularity  is  known  to  asymp- 
totically exist  along  a crack  front,  there  appears  to  be  no  additional  math- 
ematical insight  that  can  be  exploited  for  developing  a formulation  generally 
applicable  to  buried,  surface,  and  through  crack  problems.  It  cannot  be 
determined  a priori  over  what  extent  the  singularity  dominates,  and  thus, 
what  size  elements  should  be  used  at  the  front.  Hence,  strictly  speaking 
a convergence  study  should  be  a part  of  any  crack  analysis.  Whereas  for  2D 
problems,  elements  which  contain  higher  order  terms  of  the  crack  tip  expan- 
sion can  be  used  as  an  alternative  to  very  small  tip  elements  which  represent 
only  the  leading  singular  term,  this  is  not  possible  for  3D  problems,  since 
the  2D  plane  strain  crack  tip  expansion  only  holds  very  close  to  the  front. 

The  encouraging  fact  is  that  accurate  solutions  to  important  test  proble-  ; 
have  been  obtained  using  singularity  elements  with  quit*,  '^nageable  meshes, 
indicating  the  viability  of  the  approach  for  general  applicai  ons. 

The  problems  discussed  are  those  described  in  the  publications: 

Int'l.  Jour.  Fracture,  _9,  pp  340-343,  1973,  and  Nuc.  Engr.  and  Design,  26, 
pp  282-290,  1974.  In  this  work  the  assumed  displacement  method  was  used 
with  six  node  we.ige  shaped  singularity  elements  surrounding  the  crack  front. 

The  interpolation  function  had  displacements  in  planes  normal  to  the  front, 
depending  upon  the  square  root  of  distance  from  the  front  and  depending  upon 
the  angular  and  front  direction  coordinates  in  a regular  fashion.  The  dis- 
placement component  along  the  front  was  given  a form  with  non-singular 
gradients,  and  thus  the  element  was  designed  for  a combined  Mode  I-1I  situ- 
ation. As  in  standard  finite-element  formulations,  the  unknowns  were  nodal 
displacement  values.  The  singular  element  has  the  necessary  rigid  translation 
modes.  It  does  not  have  the  linear  displacement  mode  which  would  represent 
rigid  body  rotation  or  constant  strain  (such  as  for  average  thermal  expansion). 
Clearly,  the  significance  of  a linear  displacement  mode  relative  to  the 
square  root  mode  diminishes  as  the  front  is  approached  (except  for  the  trivial 
cases  where  = 0,  as  in  unconstrained  aniferm  thermal  expansion). 

Therefore,  an  accurate  convergent  solution  is  attaina".  le  with  the  e’ ement  by 
successively  decreasing  element  size.  Practically  speaking,  however,  this 
matter  of  linear  displacement  mode  can  be  a serious  constraint  on  element 
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size  and  must  be  considered  on  a case  by  case  basis. 

The  test  problems  are  the  buried  penny  crack,  the  semi-circular 
surface  crack,  the  quarter-circular  corner  crack,  and  the  compact  tension 
specimen.  The  K distributions  are  described,  and  also  the  general  features 
of  the  complete  solutions.  Comparisons  with  other  published  works  suggest 
that  the  singularity  element  solutions  are  highly  accurate.  With  the  wide 
range  of  applicability  of  the  singularity  element  approach  and  demonstrated 
accuracy,  it  is  felt  that  great  progress  has  been  made  in  our  capability  to 
analyze  3D  cracked  hardware. 


DISCUSSION 


The  discussion  following  this  presentation  was  concerned  with  the 

powers  of  r used  in  the  finite-element  solution  and  element  sizes  near  the 

crack  tip.  The  term  r^^  was  included,  but  a procedure  for  including  ra, 

where  a would  be  determined  by  the  program,  was  not  part  of  the  study. 

Tracey  reported  that  the  decreasing  element  size  gave  a worse  solution. 

1/2 

This  was  attributed  to  the  area  over  which  the  r 


term  was  active. 
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AN  EVALUATION  OF  THE 

QUADRATIC  ISOPARAMETRIC  SINGULARITY  ELEMENT 

by 

Joseph  M.  Bloom 

The  Babcock  & Wilcox  Company 
Alliance,  Ohio 


Concern  for  a reliable  technique  to  determine  stress  intensity 
factors  for  both  two-dimensional  and  three-dimensional  crack  problems  led 
the  author  to  an  evaluation  of  the  available  crack  tip  elements.  In  late 
1974  and  early  1975  the  author  became  aware  of  the  potential  of  the  qua- 
dratic isoparametric  8-node  planar  and  20-node  solid  finite-elements  as 
(1  2) 

singular  elements  ’ . The  evaluation  of  the  two-dimensional  8-node  qua- 

(3) 

dratic  isoparametric  singularity  was  made  during  1975  and  presented  at 
the  Ninth  National  Symposium  on  Fracture  Mechanics  in  Pittsburgh  on  August  8, 
1975.  An  evaluation  of  the  convergence  characteristics  of  the  three-dimen- 
sional 20-node  quadratic  isoparametric  singularity  brick  element  is  presented, 
as  well  as  a brief  study  of  the  8-node  axisymmetric  singular  element. 

The  compact  tension  specimen  was  chosen  as  the  example  problem 
to  evaluate  the  20-node  singular  element.  Several  investigators  have  ana- 
lyzed this  Standard  ASTM  E399-74  fracture  specimen^  Based  solely 

on  these  references,  it  is  not  clear  to  the  author  that  convergence  of  the 
respective  solutions  in  these  references  has  been  demonstrated.  Three  grid 
refinements  of  64,  216,  and  312  elements  were  selected  to  study  the  con- 
vergence. The  mid-surface  stress  intensity  was  determined  using  the  dis- 
placement extrapolation  method  with  the  vertical  displacements  along  the  free 
surface  of  the  crack.  Table  1 presents  the  results  for  both  the  three- 
dimensional  condition  and  the  plane  strain  condition,  in  which  displacements 
perpendicular  to  the  specimen  faces  were  set  to  zero.  The  312  element  grid 
gave  a stress  intensity  factor  to  within  1 percent  of  the  accepted  plane 
strain,  value;  while  the  three-dimensional  condition  gave  a midplane 
stress  intensity  factor  7.6  percent  greater  than  the  accepted  K value. 

/O') 

Earlier  finite  difference  results  of  Ayres ^ ' for  a center  through-thickness 
cracked  plate  suggest  an  elevation  of  stress  intensity  at  midplane  of  10 
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percent  above  1^.  More  recent  work  by  Schroedl  and  Smith''  using  the 
photoelastic  stress  freezing  technique,  gave  calculated  midplane  stress 
intensity  factors  of  from  8 to  10  percent  above  The  compact  tension 

problem  was  rerun  with  wedge  elements  surrounding  the  crack  front  and,  while 
the  extrapolation  value  of  stress  intensity  did  not  change,  the  stress 
intensity  determined  from  the  quarter-point  node  displacement  on  the  crack 
for  the  plane  strain  case  increased  from  7,008  ksi  /in.  to  7,134  ksi  /in. 
giving  a 0.9  percent  difference  compared  to  the  value  found  by  Brown 
and  Srawley.  For  the  three-dimensional  condition  with  318  elements  (wedges 
around  the  crack),  this  difference  was  found  to  be  7.8  percent  higher  than 
the  result  for  the  midplane  value  of  the  stress  intensity. 

In  addition  to  the  above  study,  the  8-node  axisymmetric  singular 
element  was  evaluated.  Two  axisymmetric  geometries  were  run:  a solid 

cylinder  with  a penny-shaped  crack  and  a cylindrical  shell  with  a circum- 
ferential flaw.  The  solid  cylinder  was  run  with  a uniform  axial  load.  The 
cylindrical  shell  was  run  with  both  a remote  uniform  load  and  a non-uniform 
load  acting  over  the  crack.  Two  crack  sizes  (a/t)  were  run  for  both  load 
conditions.  Comparison  of  stress  intensity  factors  with  accepted  solutions 
from  the  literature  showed  the  8-node  axisymmetric  singular  element  gives 
excellent  results. 

Based  on  the  present  evaluation,  it  appears  that  both  the  20-node 
singular  isoparametric  brick  element  and  the  8-node  singular  axisymmetric 
rectangular  element  are  excellent  singularity  elements  which  require  only 
a relatively  few  elements  for  obtaining  accurate  stress  intensity  factors. 

In  addition,  the  non-singular  form  of  the  element  is  currently  in  most 
f inite-element  stress  analysis  computer  programs,  thereby  eliminating  the 
ncessity  of  adding  a special  crack  tip  element  to  the  existing  computer  code. 
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TABLE  1.  COMPACT  TENSION  SPECIMEN 

20-NODE  QUADRATIC  ISOPARAMETRIC 
THREE-DIMENSIONAL  SINGULARITY  ELEMENT 


DECREES 

OF  „ 1 % * 

NO.  OF  ELEMENTS  FREEDOM  Kl  DIFFERENCE 


THREE  - DIMENSIONAL  CONDITION 


64 

1200 

7450 

+3.5 

216 

3520 

7650 

+6.3 

312 

4005 

7750 

+7.6 

1 MIDPOINT  Z - DISPLACEMENT 

2 ONE  LAYER  OF  20  - NODEO  CLEMENTS 
• COMPARISON  WITH  K2Q  = 7200  PSI  IN 


i ] 

i 

* 

i 

j 

» 

1 ^ 
i ) 


i, 


t 


i 


i 

i 


\ 

\ 


A'L: 


Ihtfjjijl^^-'^  * * Wl*,?^  * *"'’**  'v'  ” ' ' , 


20 


REFERENCES 


1.  R.  D,  Henshell  and  K.  G.  Shaw,  "Crack  Tip  Finite-Elements  are  Unnecessary", 
Research  Report,  University  of  Nottingham  (1973).  Also  in  the  Inter- 
national Journal  for  Numerical  Methods  in  Engineering,  Vol.  9,  No.  2 
(1975). 

2.  R.  S.  Barsoum,  "Application  of  Quadratic  Isoparametric  Finite-Elements 

in  Linear  Fracture  Mechanics",  International  Journal  of  Fracture,  Vol.  10 
(1974). 

3.  J.  M.  Bloom,  "An  Evaluation  of  a New  Crack  Tip  Element — The  Distorted 
8-Node  Isoparametric  Element",  International  Journal  of  Fracture,  Vol.  11 
(1975). 


fi 


4.  D.  M.  Tracey,  "Finite-Elements  for  Three-Dimensional  Elastic  Crack 
Analysis",  Nuclear  Engineering  and  Design,  Vol.  26  (1974). 

5.  T.  A.  Cruse,  "An  Improved  Boundary-Integral  Equation  Method  for  Three- 
Dimensional  Elastic  Stress  Analysis",  Computer  and  Structures,  Vol.  4, 
(1974). 

6.  R.  S.  Barsoum,  "On  the  Use  of  Isoparametric  Finite-Elements  in  Linear 
Fracture  Mechanics" , International  Journal  for  Numerical  Methods  in 
Engineering,  Vol.  10  (1976). 

7.  J.  Reyner,  "On  the  Use  of  Finite-Elements  in  Fracture  Analysis  of  Pressure 
Vessel  Components" , ASME  Paper  No.  75-PVP-20  presented  at  the  Second 
National  Congress  on  Pressure  Vessels  and  Piping,  San  Francisco,  Cali- 
fornia, June  23-27,  1975. 

8.  D.  J.  Ayres,  "A  Numerical  Procedure  for  Calculating  Stress  and  Deforma- 
tion Near  a Slit  in  a Three-Dimensional  Elastic-Plastic  Solid",  Engineer- 
ing Fracture  Mechanics,  Vol.  2 (1970). 

9.  M.  A.  Schroedl  and  C.  W.  Smith,  "Influence  of  Three-Dimensional  Effects 
on  the  Stress  Intensity  Factors  of  Compact  Tension  Specimens",  ASTM 
STP  560  (1974). 


I 

I 

I 


WP??- 


21 

DISCUSSION 


The  boundedness  of  the  strain  energy  density  was  discussed. 
Bloom  briefly  showed  the  basis  for  the  existence  of  a bounded  energy.  He 
found  that  triangular  elements  gave  better  results  than  rectangular  ele- 
ments, but  if  one  is  extrapolating,  it  does  not  make  a big  difference. 
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ON  SOME  FRACTURE  ANALYSIS  RESULTS 
AT  THE  BERKELEY  NUCLEAR  LABORATORIES 

by 

C.  H.  A.  Townley 

Central  Electric  Generating  Board 
Berkeley  Nuclear  Laboratories 
Berkeley,  Glouster shire,  England 


This  is  an  extemporaneous  address  on  some  investigations  being 
carried  out  primarily  by  T.  Hellen  and  W.  Blackburn  and  their  colleagues. 

With  the  objective  of  solving  three-dimensional  crack  problems, 
an  initial  study  was  carried  out  to  carefully  compare  results  obtained  in 
two-dimensional  problems  with  four  approaches. 

1.  Standard  analysis  with  standard  and  special  crack-tip 
finite- elements 

2.  Solving  crack  problems  with  two  slightly  different  length 
cracks 

3.  Virtual  crack  extension  approach 

A.  J- integral  approach,  modified  for  thermal  stress. 

Comparisons  between  these  methods  were  discussed.  It  appeared 
that  special  crack  tip  elements  were  necessary  in  most  cases  to  obtain  satis 
factory  solutions  with  all  of  the  approaches. 

One  three-dimensional  crack  considered  was  the  cross-corner  crack 
in  a pipe  T-connection.  This  problem  was  investigated  with  a finite-element 
program  with  substructuring,  in  which  a special  cracked  substructure  was 
developed.  Four  different  crack  sizes  were  solved. 

A second  three-dimensional  problem  studied  was  the  "pop-in"  pro- 
blem for  compact  tension  specimens  with  "hard"  loads.  This  problem  was 
analyzed  both  for  a straight-through  crack  and  for  a thumbnail  crack. 
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DISCUSSION 

There  was  a general  discussion  following  the  last  presentation. 
Discussions  of  some  topics  brought  up  are  summarized  here.  The  question  of 
the  t-'pe  of  behavior  in  the  isoparametric  element  with  nodes  at  the  quarter 
points  was  brought  up  by  J.  Bloom  with  the  specific  point  that  there  are 
papers  in  the  literature  contending  that  this  type  of  element  is  not  correct 
because  the  strain  energy  is  not  valid.  Bloom  pointed  out  that  the  extra- 
polation procedure  produces  results  similar  to  other  analyses. 

The  question  of  the  behavior  in  the  compact  tension  specimen  was 
discussed  by  T.  Cruse.  It  appears  that  some  data  show  that  results  converge 
to  the  plane  strain  value  of  K while  others  converge  to  10  percent  greater 
than  that  value.  Various  aspects  of  the  solution  technique  including  the 
type  of  finite-element,  the  variation  of  displacements  through  the  thickness 
and  the  sensitivity  of  the  solution  to  small  changes  in  the  loading  were 
discussed.  It  was  brought  out  that  the  comparisons  being  discussed  have  not 
been  defined  precisely  enough  to  determine  if  the  problems  being  solved  are 
indeed  the  same  ones.  J.  Swedlow  pointed  out  that  the  sensitivity  of  the 
compact  tension  solution  to  changes  in  the  loading  requires  that  loadings  be 
described  in  fine  detail  in  order  for  valid  comparisons  to  be  made.  T.  Cruse 
said  that  he  foun-  a variation  of  10  percent  in  the  resulting  K for  a compact 
tension  specimen  just  by  introducing  the  loads  in  a different  way. 

G.  Sih  brought  out  the  point  that  three  types  of  numericrl  solu- 
tions for  the  value  of  K approaching  a free  surface  were  presented.  In  one 
case,  it  goes  up,  in  another  it  goes  up  and  comes  down,  and  the  the  third 
case,  it  goes  down.  His  comments  were  that  different  numerical  approaches 
can  contribute  to  the  characteristics  of  the  solution  and  that  one  should 
not  discuss  5 to  10  percent  differences  until  the  techniques  are  defined. 

While  it  is  not  known  what  the  three-dimensional  solutions  should  be  it  was 
generally  agreed  that  for  two-dimensional  solutions,  correlation  with  plane 
strain  results  is  necessary,  but  not  sufficient  for  verification  of  an  analysis 
technique. 

(X 

Reasons  for  introducing  a term  of  the  type  r jn  the  solution  and 
questions  about  the  value  of  a being  different  on  the  free  surface  and  the 
interior  were  discussed.  These  discussions  as  well  as  previous  discussions, 
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pointed  out  that  the  present  state  of  knowledge  about  the  three-dimensional 
solution  is  such  that  the  characteristic  power  of  r at  the  free  surface 
is  not  known  and  while  this  question  is  being  addressed,  an  exact  elasticity 
solution  does  not  appear  to  be  on  the  immediate  horizon.  Concerning  numerical 
techniques,  observations  that  the  form  of  the  strain  energy  within  the  quarter 
point  isoparametric  elements  loses  it's  r characteristic  were  expressed. 
However,  numerical  results  pointing  to  convergence  of  numerical  solutions 
and  theoretical  bases  for  convergences  were  brought  out  to  support  the  utility 
of  the  numerical  techniques.  Thus,  while  the  three-dimensional  numerical 
techniques  are  a valuable  tool  for  the  stress  analysis,  there  are  no  exact 
three-dimensional  solutions  tc  serve  as  a basis  for  comparison.  This  pro- 
vides a fertile  ground  for  discussion  and  future  research  in  this  area. 
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STRESS  INTENSITY  ESTIMATES  FOR  THREE-DIMENSIONAL 
CRACKED  BODY  PROBLEMS  BY  THE  FROZEN  STRESS 
PHOTOELASTIC  METHOD 

by 

C.  W.  Smith 

Virginia  Polytechnic  Institute  and  State  University 
Blacksburg,  Virginia 
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Based  upon  an  idea  of  G.  R.  Irwin ^ and  analytical  studies  of 

(2) 

G.  C.  Sih  and  his  associates,  the  senior  author  and  his  associates  have 

(3-2?) 

developed,  over  a period  of  years  a computer  assisted,  frozen  stress 

photoelastic  technique  for  estimating  stress  intensity  factors  (SIFs)  in 
three-dimensional  cracked  body  problems.  Originally  developed  for  Mode  I 
problems  only,  the  method  has  been  extended  to  Mixed  Mode  problems  (i.e., 
Mode"!  I and  II)  and  studies  are  currently  being  conducted  which  are  directed 
towards  the  inclusion  of  Mode  III  effects  in  the  measurements. 

For  Mode  1 SIF  estimates,  the  analysis  involves  expressing  the 
maximum  shearing  stress  in  a plane  mutually  perpendicular  to  the  flaw  border 
and  the  flaw  surfaces  in  the  form 
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along  a line  normal  to  the  flaw  surfaces  and  passing  through  the  crack  tip. 
By  truncating  the  Taylor  Series  to  its  first  term  only,  and  combining  the 
result  with  the  stress-optic  law  in  two  dimensions 
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where  n’  is  the  photoelastic  fringe  order,  f is  the  material  fringe  value 
and  t'  the  slice  thickness,  an  expression  of  t'i~  form 

,1/2 


K 


AL. 

o(jia) 


1/2 


o(ira) 


1/2 


+ 


V8>- 


r_ 

a 


1/2 


(3) 


* Where  r is  distance  from  crack  tip. 
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can  be  obtained,  where  K,  - x (8irr)'w  ^ which  shows 

Ap  max 
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to 


be  a linear  relation.  By  obtaining  photoelastic  data  and  using  least  squares 

V r“  *«  I/2 


to  fit  a line  to  this  data,  the  graph  of 
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vs 


r_ 

a 


can  be  extrapo- 


lated to  the  origin  to  obtain  K,.  By  working  directly  with  x , stress 

I max 

separation  methods  are  avoided. 

Following  the  ideas  of  Reference^,  the  authors  have  developed  a 
similar  approach  for  estimating  both  K^.  and  values  in  a Mixed  Mode 
problem. 


The  experimental  data  are  obtained  by  constructing  a scale  model  of 
the  body  from  fringe  free  transparent  photoelastic  material,  inserting  a 
"natural"  or  an  artificial  flaw,  and  heating  the  structures  to  the  critical 
temperature  of  the  material.  At  this  temperature,  the  material  is  fully 
elastic,  and  it  is  loaded  to  a desired  level.  The  model  is  then  cooled  under 
load,  "fixing"  both  the  photoeiastic  fringes  and  the  deformations  so  that 
both  are  retained  in  the  material  even  after  unloading  and  slicing  of  the 
model.  Thin  slices  mutually  perpendicular  to  both  the  flaw  border  and  the 
flaw  surfaces  are  then  removed  at  designated  locations  along  the  flaw  border, 
immersed  in  a liquid  of  the  same  refractive  index  as  the  model  material  and 
analyzed  in  a crossed  circular  polariscope  via  the  Tardy  Method.  Results  are 
fed  into  a least  squares  digital  computer  program  for  estimating  the  S1F  for 
each  slice. 

Figure  1 shows  typical  fringe  patterns  for  both  Mcdc  I and  Mixed 
Mode  loadings  and  Figures  2 and  3 show  typical  results  for  two  three-dimen- 
sional problems. 

These  results  show  that,  even  for  rather  complex  three-dimensional 
problems,  the  linear  zone  can  be  located  experimentally  and  results  such  as 
those  shown  here  can  be  replicated  to  within  + 5 to  7 percent. 

Since  Poisson's  Ratio  - 0.5  for  all  stress  freezing  materials  above 
critical  temperature,  it  ir<  necessary  to  include  a correction  when  it  is 
desired  to  apply  results  to  a material  for  which  v = 0.3.  This  correction 
ranges  from  - +5  percent  for  surface  flaws  to  as  much  as  12  percent  for 
through  cracks  in  plates  of  finite  thickness.  Studies  are  in  progress  to 
accurately  quantify  this  effect. 
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To  date  a fairly  broad  variety  of  technologically  important  three- 
dimensional  problems  have  been  studied.  It  is  felt  that  the  method  offers 
a reasonable  independent  experimental  check  on  finite-element  results  pro- 
vided experimental  error  of  the  order  of  5 percent  can  be  tolerated. 
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DISCUSSION 

Swedlow  raised  an  important  point  in  that  most  photoelastic  anal- 
ysis work  is  focused  on  obtaining  the  stress- intensity  factur,  but,  other 
information  can  be  obtained  of  value  for  checking  a finite-element  compu- 
tation. Smith  concurred,  but  pointed  out  that,  while  numbers  such  as  the 
COD's  can  be  obtained  from  a photoelasticity  solution,  they  are  not  un- 
equivocal. A correction  must  be  introduced  which  depends  on  the  modulus 
and  Poisson's  ratio  of  the  material.  While  this  is  straightforward  in  two- 
dimensional  problems,  it  is  not  in  three-dimensions. 
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OBSERVATIONS  OF  CRACK  TIP  PROCESSES 

by 

G.  T.  Hahn 
BATTELLE 

Columbus  Laboratories 
Columbus , Ohio 


Several  metallographic  methods  can  be  used  to  study  three-dimen- 
sional crack  problems  in  addition  to  the  photoelastic  method  discussed 
here  by  Professor  C.  W.  Smith.  The  Fe-3Si  alloy-dislocation  etching  tech- 
nique^ has  been  used  to  good  advantage  to  reveal  the  plastic  zone  of 
cracks.  Figure  1 shows  examples  of  the  plastic  zone  produced  by  a fatigue 

crack  both  on  the  surface  and  in  the  interior  on  the  midsection  of  a compact 

(2) 

specimen  . This  etching  technique  is  specific  to  Fe-3Si  steel,  a material 
with  a stress  strain  curve  similar  to  medium  strength  constructional  steels. 
It  has  been  used  to  delineate  crack  tip  plastic  zone  produced  by  monotonic 
loading^,  cyclic  loading  stable  crack  growth^,  and  fast  fracture^ 

Three-dimensional  effects  have  an  important  role  in  the  evaluation 
of  the  material  properties  governing  fast  fracture  and  crack  arrest.  Figure 
2 shows  a.i  example  of  the  profile  of  an  arrested  fast  fracture  in  the  A533B 
steel  test  section  of  a duplex  DCB  test  piece v . Figure  2a  (about  13  mm 
from  the  crack  tip)  shows  microscopic  branch  nuclei,  which  are  believed  to 
be  responsible  for  the  macroscopic  branching  event  observed  in  a comparison 
specimen  shown  in  Figure  3.  The  branch  nucleation  event  is  a poorly  under- 
stood three-dimensional  phenomena  which  complicates  the  measurement  and 
interpretation  of  propagation  and  arrest  events.  Deep  side  grooves,  which 
can  be  used  to  restrict  branching  events  of  this  kind  also  introduce  a three- 
dimensional  component  to  the  stress  field  of  the  crack.  Figures  2a  and  2b 
show  ligaments,  which  are  another  three-dimensional  feature  of  predominatly 
cleavage  fractures  There  is  evidence  that  the  major  part  of  the  tough- 

ness displayed  by  steels  tested  below  the  transition  temperature  is  derived 
from  the  ductile  fracture  of  the  ligaments  which  contributed  a small  fraction 
of  the  total  fracture  surface^. 
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Figure  4 shows  the  arrested  crack  front  which  was  revealed  by 
heat  tinting  the  crack  surface  (dark  portion)  before  breaking  the  specimen 
apart  at  -196°  C (light  portion).  The  light  line  visible  behind  the  crack 
front  is  a long  unbroken  ligament  (reverse  tunnel)  which  exerted  a drag  on 

/ON 

the  crack  front.  ' Analyses  that  can  relate  such  perturbations  of  the 
crack  front  to  local  variations  in  the  fracture  resistance,  would  be  useful. 
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FIGURE  4.  PHOTOGRAPH  OF  AN  ARRESTED  CRACK  FRONT.  The  crack  front, 
which  propagated  from  left  to  right,  was  produced  by  wedge 
loading  a A533B  steel  DCB-test  piece  (with  deep  sidegrooves) 
at  -78  C.  The  front  was  delineated  by  heat  tinting  the 
surface  of  the  arrested  crack  (dark  area)  and  then  breaking 
the  specimen  open  (light  fracture). 
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DISCUSSION 

Cruse  initiated  an  extensive  discussion  regarding  tunneling 
(where  the  crack  advances  faster  in  the  inerior  than  on  the  surfaces)  and 
reverse  tunneling  as  a three-dimensional  phenomenon.  No  conclusions  were 
reached  as  to  how  it  might  be  treated,  however. 


USE  OF  CYCLIC  GROWTH  TESTS  TO  INFER  STRESS  INTENSITY 

by 

J.  E.  Collipriest 

Rockwell  International,  Space  Division 
Downey,  California 


The  technique  of  inferring  stress  intensity  values  from  experi- 
mental crack  growth  tests  aM  a knowledge  of  a crack  growth-stress  intensity 
relationship  is  increasing  in  popularity  as  a tool  for  treating  complex 
geometry  problem.1.  Some  caution  is  required  to  insure  that  the  inferred 
stress  intensity  is  correctly  depicted  and  not  confused  by  material  and/or 
experimental  variables. 

For  illustration,  the  preliminary  assessment  of  some  test  data  is 
discussed.  This  data  came  from  constant  amplitude  crack  growth  tests  of 
center  cracked  tension  geometry  specimens  with  various  sized  through  holes 
centrally  located.  These  specific  through-crack  tests  were  conducted  to 
verify  the  experimental  technique  and  material  behavior  prior  to  conducting 
tests  of  holes  with  corner  cracks. 

The  initial  data  assessment  showed  three  anomalies  when  inferred 
stress  intensity  values  were  compared  to  theoretical  values  for  the  CCT  geom- 
etry with  the  Bowie  treatment  for  near  hole  effects  and  the  Fedderson  secant 
correction  for  finitie  width  effects.  For  crack  lengths  up  to  twice  the 
hole  radius,  there  appeared  to  be  a stress  level  effect  wherein  the  lower 
the  cyclic  stress  the  higher  the  implied  stress  intensity  factor.  For  long 
cracks,  tne  stress  intensity  factors  started  to  increase  more  sharply  than 
the  secant  width  correction  would  provide.  And  finally,  there  were  several 
indications  of  a random  shift  of  implied  stress  intensity  when  higher  growth 
rates  were  experienced. 

The  short  crack  length  anomaly  is  shown  to  be  the  result  of  a 
subtle  mis-representation  of  cyclic  growth  rate-vs -delta  K data  for  growth 
rates  below  10  micro-inches  per  cycle. 

The  long  crack  anomaly  can  be  numerically  accommodated  by  arbi- 
trarily adjusting  the  secant  width  correction  by  10  percent.  The  only  basis 
suggested  for  this  arbitrary  adjustment  is  that  it  reasonably  fits  the  data 
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for  a fairly  large  variety  of  specimen  widths  and  crack  lengths  greater 
than  30  percent  of  specimen  width.  Additionally,  it  was  noted  that  a similar 
adjustment  to  the  secant  width  correction  was  required  to  match  elastic  com- 
pliance data  for  CCT  specimens  with  long  crack  lengths. 

The  third  anomaly  is  attributed  to  the  variety  of  three-dimensional 
geometries  that  through  cracks  may  assume  at  high  stress  intensities  and 
high  growth  rates.  At  the  transition  from  flat  to  slant  crack  geometry  two 
basic  forms  may  occur:  full  slant  or  cup-cone.  Since  the  center  crack  geom- 
etry has  two  crack  fronts  a large  variety  of  combinations  will  occur  such  as 
matching  or  opposed  slant,  matching  or  opposed  cup-cone  and  combinations 
where  one  crack  is  slant  and  the  other  is  cup-cone. 

Substantial  variations  of  out  of  plane  displacements  result  and 
the  differing  resultant  mixed  mode  stress  intensities  are  reflected  in 
differing  crack  growth  rates.  All  of  which  can  be  avoided  by  only  conducting 
tests  at  low  stress  intensity  values  so  that  cracks  are  only  grown  in  a flat 
plane  normal  to  the  principal  applied  stress. 

One  further  complicating  factor  was  illustrated  by  crack  growth 
data  for  one  aluminum  alloy,  single  product  form,  single  thickness.  Signif- 
icant variation  of  growth  rates  were  observed  between  material  from  different 
manufacturers.  A actor  of  two  variation  occurred  at  10  ksi  /in  and  5 micro- 
inches/cycle and  more  than  an  order  of  magnitude  difference  occurred  above 
20  ksi  /in. 

It  was  concluded  that  the  technique  of  inferring  stress  intensity 
from  cyclic  growth  rate  tests  was  viable.  However,  considerable  care  is 
required  in  test  conduct  and  reconciliation  of  results. 
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DISCUSSION 

Cruse  was  disturbed  by  the  fact  that  obvious  large  cyclic  plas- 
ticity and  the  possibility  of  mixed  modes  were  present  in  the  experiments, 
yet  comparisons  were  made  from  specimen  to  specimen  where  significant  dif- 
ferences in  these  conditions  were  ignored.  Specifically,  the  fatigue  crack 
growth  data  involved  plasticity  and  mixed  mode  mechanisms  that  are  necessarily 
thickness  dependent.  Correlation  can  be  expected  only  where  the  plate  thick- 
ness and  mode  of  crack  growth  are  identical.  On  this  basis,  Cruse  suggested 
that  plasticity  could  play  a large  role  in  the  work  and  that  the  correlation 
problems  evidenced  by  Collipriest  might  be  explained  in  this  way. 
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SUMMARY  OF  COMMENTS  FOR  WORKSHOP  ON  THREE-DIMENSIONAL 
FRACTURE  ANALYSIS 

by 

J.  L.  Swedlow 

Carnegie-Mellon  University 
Pittsburg,  Pennsylvania 


It  is  evident  that  analytical  (computational)  procedures  are  now 
in  hand  to  attack  a useful  range  of  three-dimensional  elastic  problems  in- 
volving cracks.  It  is  my  sense,  however,  that  individual  cases  selected  for 
study  are  frequently  chosen  more  on  the  basis  of  the  simplicity  of  their 
geometrical  and  loading  configuration  than  the  utility  of  the  end  result. 

That  is,  the  issues  encountered  in  laboratory  (data-collecting)  and  service 
(data-applying)  situations  should  be  more  influential  when  problems  are 
posed  for  analysis.  To  exemplify  this  view,  one  situation  of  some  pertinence 
is  reviewed. 

(1  2) 

Some  experimental  information  ' has  been  developed  which  sug- 
gests that,  at  least  In  a compact  specimen,  crack-tip  K values  are  not  what 

(3) 

would  be  computed  using  standard  formulae,  e.g.,  . While  the  nature  and 

scope  of  this  work  does  not  yet  give  a full  picture  of  behavior  germane  to 
the  issue,  two  notions  emerge  from  even  these  limited  data: 

(1)  Formulae  for  K pertinent  to  the  compact  specimen  are 
especially  sensitive  to  details  of  load  arrangement 

(as  opposed  to  overall  magnitude).  This  point  is  noted 
elsewhere ^ and  may  carry  over  to  "three-dimensional" 
configurations  currently  of  popular  interest,  e.g.,  the 
corner  crack  embedded  in  a pressurized  vessel. 

(2)  Configuration  of  the  crack  itself  may  have  a dispropor- 
tionate influence  on  local  K values.  An  early  suggestion 
to  this  effect ^ concerns  a number  of  geometries  other 

. ,-an  a standard  - or  nearly  so  - test  specimen.  Some 

dramatic  results  obtained  in  the  convenient  cases  were 
(5) 

treated  in  . 
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The  relevance  of  these  notions  to  three-dimensional  fracture  analysis  is  to 

be  noted  in  two  quite  different  respects.  Both  show  that  accurate  K values 

are  not  of  necessity  derived  from  planar  analysis,  either  in  some  global  or 

average  sense,  or  along  the  crack  front.  Were  such  information  in  hand, 

however,  interpretation  of  some  test  results  would  be  clarified.  As  noted 

in^,  for  example,  the  simple  load-carrying  capacity  of  a specimen  derived 

(3  4) 

from  analyses  of  the  form  summarize!  in  sources  such  as  * may  in  some 
instances  become  misleading.  The  inference  drawn  in^  is  that  interpreta- 
tion of  fracture  test  data  may  be  nonconservative.  Of  perhaps  greater 
importance  is  the  collection  and  use  of  fatigue  test  data.  Under  cyclic 
loading  of  a CS  specimen,  the  crack  grows  and  the  contours  of  its  front  may 
alter;  no  cognizance  is  taken  of  this  matter  in  computing  a planar  AK,  which 
is  the  sole  mechanism  for  transferring  data  from  the  laboratory  to  design  or 
service  situations. 

The  relevance  of  better  K values  along  the  cxack  front  is  also 
seen  in  three-dimensional  fracture  analysis,  the  subject  of  this  workshop. 

At  present,  there  is  less  than  a clear  view  of  what  problems  need  attention. 
Manifestly,  numerical  procedures  have  been  developed  to  the  point  where  a 
wide  range  cf  problems  may  be  solved,  but  the  setting  of  actual  cases  to  be 
treated  seems  to  be  guided  only  incidently  by  physical  observation.  While 
the  results  reported  here  represent  major  Investments  of  enevgy  and  intellect, 
there  remains  need  for  close  attention  to  actual  behavior  of  a series  of 
specific  shape/material/loading  systems  so  that  the  computational  power  now 
available  is  advantageously  used. 
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DISCUSSION 


The  data  of  Neale,  introduced  by  Swedlow,  related  the  ratio  of 
the  local  value  to  the  nominal  K as  a function  of  the  crack  front  curvature. 
Hahn  noted  that  these  data  appear  to  conflict  with  the  Barsoum-Clausing 
results  in  which  positive  tunneling  increases  the  K level.  Neale's  data 
gives  a linear  decrease  with  crack  fron  curvature.  Townley  suggested 
that  viewing  the  appearance  of  a thumbnail  as  a stability  effect  Piay  be 
helpful. 
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CRACK  TIP  FIELDS  IN  STEADY  CRACK  GROWTH 
WITH  STRAIN  HARDENING 

by 

J.  W.  Hutchinson 

Harvard  University 
Cambridge,  Massachusetts 


Singular  stress  and  strain  fields  are  found  at  the  tip  of  a crack 
growing  steadily  and  quasi-statically  into  an  elastic-plastic  strain  hard- 
ening material.  'The  material  is  characterized  by  flow  theory  together 
with  a bilinear  effective  stress-strain  curve.  Anti-plane  shear,  plane 
stress  and  plane  strain  are  each  considered.  Numerical  results  are  obtained 
for  the  stress  and  strain  fields,  the  order  of  the  singularity  and  the  near 
tip  regions  of  plastic  loading  and  elastic  unloading. 
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SOME  PROPERTIES  OF  FINITE-ELEMENT  APPROXIMATIONS  OF 

ELLIPTIC  PROBLEMS  ON  DOMAINS  WITH  CRACKS  AND  CORNERS  ! 

f 

by  | 

i 

J.  T.  Oden 

The  University  of  Texas 
Austin,  Texas 

t 

A survey  of  some  of  the  mathematical  foundations  of  elliptic 
variational  theory  for  two-dimensional  domains  with  corners,  including 
cracks,  is  described.  The  basic  mission  here  is  to  describe  the  mathe- 
matical framework  in  which  analyses  of  finite-element  methods  for  problems 
in  fracture  mechanics  must  be  studied.  In  particular,  the  theory  of 
weighted  Sobolev  spaces  is  described  together  with  appropriate  imbedding 
theorems  and  an  existence  theorem  for  variational  boundary-value  problems 
set  in  these  spaces.  It  is  shown  how  these  results  help  form  the  basis 
for  obtaining  error  estimates  for  restricted  classes  of  finite-element 
methods.  A priori  error  estimates  in  energy  norms  in  an  L^-norm  and  in  the 
L^-  norm  are  described.  The  study  summarizes  the  principal  results  of  , 

Babuska  and  some  very  recent  findings  of  Schatz  and  Whalbin.  ! 
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NEAR  FIELD  BEHAVIOR  AND  CRACK  GROWTH 


by 

George  C.  Sih 

Institute  of  Fracture  and  Solid  Mechanics 
Lehigh  University 
Bethlehem,  Pennsylvania 


A knowledge  of  the  local  three-dimensional  stress  field  for  a 
crack  with  an  arbitrarily  curved  front  is  necessary  for  investigating  the 
crack  growth  behavior.  Such  a stress  field  is  referred  to  a system  of 
local  spherical  coordinates  (r,  6,  <j>)  and  the  result  can  be  reduced  to  a 
surprisingly  simple  form  when  the  appropriate  choice  of  coordinates  (r , 

0,  10)  are  used. 

Crack  growth  directions  for  various  positions  along  the  crack 
front  can  be  determined  from  the  strain  energy  density  fracture  criterion. 
The  mixed  mode  loading  of  a flat  elliptical  crack  serves  as  one  of  the 
examples.  The  development  of  thumbnail  cracks  in  a thick  plate  can  also 
be  predicted.  A future  application  of  this  result  to  the  fracture  thickness 
problem  is  the  interaction  of  a curved  (thumbnail)  crack  front  with  the 
highly  distorted  or  yielded  zones  (shear  lips)  near  the  plate  surfaces.  The 
analysis  requires  crack  growth  under  mixed  mode  loading  as  the  shear  lips 
are  developed  on  planes  inclined  to  the  plate  surfaces. 
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STRESS  INTENSITY  FACTOR  MEASUREMENTS 
FOR  CORNER  CRACKED  HOLES 

by 

A.  F.  Grandt,  Jr. 

Air  Force  Materials  Laboratory 
Wright-Patterson  A.  F.  B. , Ohio 


The  fatigue  crack  growth  rate  method  was  used  to  measure  stress 
intensity  factors  at  various  points  along  the  border  of  a corner  crack  lo- 
cated at  the  edge  of  hole  in  a plate  loaded  in  remote  tension.  Cyclic 
extension  of  the  corner  flaws  was  recorded  by  time  lapse  photography  in 
plexiglass  test  specimens.  Since  the  specimens  were  transparent,  it  was 
possible  to  photograph  full  plan  views  of  the  part  through  crack  and,  thus, 
determine  the  variation  in  fatigue  crack  growth  rate  around  the  flaw  peri- 
meter. The  measured  fatigue  crack  growth  rates  were  then  used  with  the 
Paris  relation  between  fatigue  c/ack  growth  rate  and  range  in  stress  inten- 
sity factor  (previously  established  for  the  test  material)  to  compute  the 
cyclic  range  in  stress  intensity.  The  stress  intensity  factors  at  the 
points  where  the  crack  intersects  the  edge  of  the  hole  and  the  front  surface 
of  the  plate  are  then  compared  with  various  analytical  predictions  avail- 
able in  the  literature.  Computations  for  other  intermediate  points  along 
the  crack  boundary  are  in  progress  and  will  be  discussed  as  available  at 
the  time  of  the  meeting. 
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DEVELOPMENT  OF  PROCEDURES  FOR  ANALYZING  STRESSES 
IN  CRACKED  BODIES  OF  VARIOUS  SHAPES 

by 

J.  C.  Bell 
BATTELLE 

Columbus  Laboratories 
Columbus,  Ohio 


During  the  last  few  years,  an  effort  has  been  made  at  Battelle's 
Columbus  Laboratories  to  develop  procedures  for  analyzing  crack  stresses 
in  bodies  involving  progressively  more  detail  of  body  shape.  The  progress 
that  has  been  made  is  built  on  two  basic  forms  of  stress  analysis  and  on 
a procedure  for  merging  analyses  of  these  forms  in  unified  calculations. 
Recent  efforts  have  centered  around  choosing  details  in  the  merger  process 
so  that  dependable  analyses  will  be  obtained.  Various  aspects  of  the  work 
have  been  .ponsored  by  Battelle,  by  the  Air  Force  Flight  Dynamics  Laboratory, 
and  by  the  NASA,  Lewis  Research  Center. 


BASIC  ANALYSES 


The  basic  analyses  designed  for  treating  cracks  is  one  for  stresses 
and  displacements  around  an  arbitrarily  loaded  circular  crack  in  an  infinite 
body.  Normal  loads  as  well  as  radial  and  circumferential  tangential  loads 
are  considered  in  terms  of  load  coefficients  which  are  arbitrarily  definable 
so  that  they  can  eventually  accommodate  body-surface  effects.  The  analysis 
for  the  circular  crack  alone  leads  to  stress  and  displacement  solutions 
expressed  in  terms  of  integrals  of  products  of  Bessel  functions.  These 
integrals  appear  formidable  at  first,  but  they  can  be  reformulated  so  as  to 
become  readily  computable  in  modern  computers.  (Thus,  a set  of  integrals 


for  a fixed  radius  and  a fixed  distance  from  the  crack  plane,  but  with 
about  1000  combinations  of  defining  indices,  can  be  comput-'i  in  about  one 
second  by  the  CDC  Cyber  73.)  The  stress-intensity  factors  of  any  mode  along 
the  crack  front,  which  are  useful  in  predicting  whether  further  cracking 


will  occur,  are  related  quite  simply  to  the  crack-load  coefficients.  Two 
papers  have  been  written  on  this  subject,  one  deriving  the  theory  (this 
has  been  submitted  for  publication) , and  one  summarizing  and  illustrating 
it  for  simple  cases  (this  has  been  accepted  for  publication  in  the  Journal 
of  Structural  Division  of  the  ASCE) . 

The  second  basic  analysis  determines  stresses  and  displacements 
in  a half-space  subjected  to  certain  elemental  loads,  both  normal  and 
tangential,  on  the  surface.  The  elemental  loads  are  quasi-pyramidal  and 
are  applied  in  overlapping  fashion  so  that  the  overall  load  on  the  surface 
is  continuous  and  varies  linearly  in  both  directions  in  each  rectangle  of 
a set  covering  the  entire  loaded  region.  The  continuity  of  the  applied 
load  pattern  avoids  the  spurious  discontinuities  of  stress  that  arise  in 
forms  of  analysis  which  use  discontinuous  (step  function)  surface  loads. 

The  expressions  for  the  stresses  and  displacements  in  the  body  are  simple, 
and  the  summation  processes  that  are  involved  are  readily  adaptable  for 
computers.  The  implied  continuity  of  stresses  on  the  surface  has  a potential 
advantage  in  that  it  allows  much  freedom  in  pointwise  fitting  boundary  con- 
ditions on  the  surface  and  meaningful  evaluation  of  boundary  condition 
residual  errors  at  nonfitted  points. 


MERGING  OF  CONTRIBUTORY  ANALYSES 


For  analysis  of  stress  around  a circular  or  part-circular  crack 
in  a finite  body,  an  analysis  of  the  first  basic  kind  is  merged  with  one  or 
more  analyses  of  the  second  kind.  This  is  accomplished  numerically  using 
a computer  program  constructed  rrom  the  formulas  of  both  kinds  of  contrib- 
utory analyses,  and  the  coordinate  transformations  needed  to  relate  them 
to  each  other.  In  broad  terms  the  merger  is  effected  by  using  the  boundary- 
point-least-squares  technique  to  find  a set  of  crack  and  surface  load  con- 
stants which  lead  to  approximate  satisfaction  of  a chosen  set  of  boundary 
conditions.  The  fitting  is  accomplished  by  a single,  least-squares  cal- 
culation, not  requiring  the  iteration  needed  by  other  investigators  who 
have  used  the  alternating  technique  for  corresponding  mergers  of  analyses. 
Once  the  proper  load  constants  have  been  found,  the  program  can  be  used  to 


evaluate  not  only  stress  intensity  factors,  but  also  stress  and  displacement 
components  anywhere  in  the  body. 

The  finding  of  load  constants  is  straightforward  in  principle, 
but  there  are  many  choices  to  be  made  in  arranging  the  calculations.  These 
include  how  long  the  series  expansions  of  crack  functions  should  be  (this 
is  decided  in  advance),  how  detailed  the  lattices  used  for  surface-load 
decomposition  should  be,  where  boundary  conditions  should  be  imposed,  and 
whether  equilibrium  should  be  enforced  among  the  surface  loads  used  to 
"free"  the  surfaces.  Trial  calculations  ha^e  shown  that  all  of  these  choices 
significantly  affect  the  computed  stress  intensity  factors,  and  since  it  is 
not  feasible  to  pursue  all  of  these  matters  into  endless  detail,  there  are 
many  questions  to  be  answered  in  considering  whether  the  results  of  a given 
calculation  may  be  regarded  as  being  dependable. 

Much  of  the  ambiguity  of  designing  calculations  for  specific 
bodies  and  load  systems  can  be  removed  by  follow' ag  standardized  design 
procedures,  but  the  question  still  remains  whether  the  chosen  procedures 
lead  to  calculations  in  accord  with  the  actual  mechanical  behavior  of  real 
bodies.  In  order  to  gain  objectivity  in  this  matter,  it  was  decided  to 
analyze  cases  for  which  experimental  stress  intensities  were  available  and 
for  which  the  choice  of  calculation  design  involves  fairly  stern  consequences. 
The  experimental  work  used  as  a reference  showed  stress  intensities  along 
the  front  of  a circular  surface  crack  in  slabs  of  varying  thickness.*  One 
particularly  challenging  case  had  a crack  four  times  as  long  as  its  depth, 
with  the  crack  penetrating  to  a depth  0.85  times  the  slab  thickness. 

The  testing  of  calculation  design  procedures  was  begun  by  anal- 
yzing cases  of  lesser  severity  (that  is,  crack  depth  to  slab  thickness  ratio 
less  than  0.85),  but  procedures  apparently  successful  for  those  cases  were 
soon  found  to  give  highly  uncertain  results  for  the  most  severe  case.  Through 
close  examination  of  results  from  many  calculations,  however,  several  pro- 
cedures have  been  chosen  which  have  removed  much  of  the  uncertainty  in  the 
computed  stress  intensity  factors.  The  basis  for  most  of  these  procedures 
can  be  understood  in  terms  of  how  individual  crack  functions  vary  on  the 


* This  experimental  work  was  performed  by  C.  W.  Smith  of  the  Virginia 
Polytechnic  Institute  under  contract  to  the  Air  Force  Flight  Dynamics 
Laboratory. 
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crack  and  in  the  slab  faces.  One  important  principle  seems  to  be  that 
crack  functions  should  be  used  to  fit  only  crack-face  conditions  and  sur- 
face loads  to  fit  only  slab-face  conditions,  except  for  the  bonafide  inter- 
actions between  the  different  load  systems.  Achieving  the  right  degree 
of  independence  of  effects  from  the  different  load  systems  while  performing 
a united  least-square  fit  is,  therefore,  one  of  the  most  important  features 
of  a good  calculation  design  procedure.  This  principle  is  especially 
important  in  the  selection  of  boundary  conditions  to  be  applied.  It  is 
possible  also  to  do  much  tailoring  of  the  crack-function  series  in  advance. 
The  design  of  surface  lattices  is  still  subject  to  ambiguities,  but  here 
also,  matters  such  as  the  interplay  between  lattice  design  and  overall 
equilibrium  have  been  considered. 

An  interesting  feature  of  some  of  the  calculated  results  is  that 
they  show  maximum  stress  intensity  factors  at  points  well  removed  from  the 
root  of  the  crack.  This  is  similar  to  the  referenced  experimental  results 
and  may  help  to  explain  the  terdency  of  circular  cracks  to  ellipticize  as 
they  grow. 


POSSIBLE  APPLICATIONS  OF  THIS  ANALYSTS 


The  computer  program  that  has  been  used  in  performing  the  calcu- 
lations has  many  features  which  have  not  yet  been  applied.  It  allows  for  the 
presence  of  six  faces,  so  that  interactions  may  be  analyzed  between  a crack 
and  the  sides  and  ends  of  the  body  nay  be  considered,  as  well  as  interactions 
from  the  top  and  bottom  faces.  Tilting  of  the  crack  is  allowable,  and  it 
may  penetrate  more  than  one  face  of  the  body.  For  consideration  of  a body 
bounded  in  three-dimensions,  the  program  allows  for  imposition  of  equili- 
brium of  overall  forces  and  moments.  Since  consideration  of  several  faces 
may  be  expected  to  make  calculations  cumbersome,  the  program  also  includes 
means  for  condensing  the  calculations  in  view  of  symmetries  which  nay  exist 
in  the  particular  body  being  analyzed. 

The  computer  program  is  complete  also  in  the  sense  that  it  contains 
formulas  for  computing  stresses  and  displacements  throughout  a body  for 
which  the  load  constants  have  been  found.  This  makes  it  possible  to  study 
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the  variation  of  stresses  in  the  vicinity  of  the  crack,  a possibility  that 
may  help  in  understanding  other  forms  of  stress  analysis  there.  Computations 
of  displacements  also  offer  an  alternative  approach  to  finding  experimental 
checks  of  the  analysis.  One  region  in  which  detailed  study  could  be  of 
particular  interest  is  around  the  end  of  the  crack.  Some  calculations  have 
shown  stress  intensities  there  near  zero,  but  further  refinements  there 
seem  to  be  needed. 

It  was  planned  originally  that  the  analysis  of  stresses  from  loads 
on  a half-space  might  be  applied  to  a considerable  variety  of  bodies  with 
surfaces  approximatable  by  sections  of  planes.  This  could  include  shapes 
obtainable  by  joining  plates  or  by  extrusion  and  could  include  curved  junc- 
tions by  using  several  planar  strips  successively  inclined  at  relatively 
small  angles,  though  some  supporting  analysis  would  be  desirable  to  smooth 
the  joining  of  analyses.  In  another  vein,  the  basic  analyses  mentioned 
here  can  be,  and  to  some  extent  have  been,  used  as  a basis  for  constructing 
a super  finite-element  including  a crack.  The  analytic  nature  of  the  solutions 
presuming  only  crack  and  surface  load  functions  lend  themselves  to  such  a 
development. 


DISCUSSION 


Fred  Smith  commented  that  the  stress  on  the  crack  circle  exterior 
to  the  surface  must  be  such  that  the  crack  does  not  close. 

C.  W.  Smith  noted  that  in  his  experimental  work,  slices  taken 
on  opposite  sides  of  the  crack  center  line  give  the  high/low  effects  at 
successive  points  in  his  K vs  angle  curves. 

R.  Shaw  noted  that  "skinny"  ellipse  results  show  an  -U  crease  in 
the  "hump"  over  the  results  shown  here. 


BOUNDARY- INTEGRAL  EQUATION 
ANALYSIS  OF  SURFACE  CRACKS* 

by 

Thomas  A.  Cruse 

Pratt  & Whitney  Aircraft 
United  Technologies  Corporation 
East  Hartford,  Connecticut 


Introduction 


Elliptical  surface  and  corner  cracks  growing  in  mode  I symmetric 
loading  are  a major  low  cycle  fatigue  (LCF)  design  problem  for  gas  turbine 
engine  structures.  During  the  majority  of  the  LCF  life  of  such  cracked 
components  these  nearly-elliptical  cracks  maintain  modest  aspect  ratios 
(semi-axis  ratio  <3)  aid  see  no  significant  "back-face",  finite  geometry 
effects.  In  order  to  predict  the  growth  of  such  cracks  it  is  necessary  to 
predict  the  mode  I stress  intensity  factor  along  the  entire  crack  perimeter, 
denoted  K(S).  The  numerical  results  presented  herein  reflect  a systematic 
development  of  a modeling  strategy  for  such  crack  problems.  Numerical 
modeling  was  accomplished  using  the  BIE  three-dimensional  stress  analysis 
computer  program  described  in  Reference  1. 

Modeling  Strategy 

The  use  of  a general  stress  analysis  method  for  fracture  mechanics 
problems  is  not  as  accurate  as  the  use  of  analysis  methods  with  "crack  tip" 
elements,  or  other  singular  techniques.  In  three-dimensional  fracture 
mechanics  problems  in  LCF,  two  measures  of  accuracy  need  be  considered.  First 
is  the  absolute  magnitude  of  the  stress  intensity  factor  at  some  crack  front 
location;  second  is  the  distribution  of  K(S)  along  the  crack  front. 


* Research  reviewed  in  this  presentation  is  being  sponsored  by  the  Air 
Force  Office  of  Scientific  Research  (AFSC) United  States  Air  Force, 
under  Contract  F44620-74-C-0060  with  Mr.  W.  J.  Walker  (NA)  as  Air  Force 
Program  Manager. 
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BIE-generated  numerical  results  obtained  for  elliptical  buried 
cracks  have  been  reported  for  a series  of  aspect  ratios  (a/b  £ 4)  and  are 
given  in  Reference  2.  The  BIE  element  map  used  in  these  studies  is  shown 
in  Figure  1.  Values  of  stress  intensity  factor  at  increments  of  15  degrees 
elliptical  angle  are  estimated  from  the  crack  opening  displacements  of  the 
row  of  nodes  nearest  the  crack  tip.  These  values  of  K(S)  have  been  found  to 
be  in  error  by  about  8-12  percent,  as  reported  in  Reference  2.  However,  by 
dividing  the  K(S)  values  for  the  ellipse  models  by  the  values  for  the  circle, 
this  systematic  modeling  error  was  essentially  canceled.  Figure  2 compares 
the  resulting  values  of  normalized  K(S) , i.e.,  the  K-distribution,  to  the 
exact  results  for  the  buried  crack. 

The  excellent  agreement  achieved  for  buried  cracks  using  this  pro- 
cedure has  been  used  for  fracture  mechanics  results  in  surface  and  corner 
crack  problems.  For  example,  a 2:1  corner  crack  K-distribution  is  obtained 
by  dividing  the  crack  opening  displacements  for  the  corner  crack  by  the 
values  for  the  buried  crack  of  the  same  aspect  ratio  (for  data  reduction 
simplicity).  The  actual  values  of  K(S)  are  obtained  by  multiplying  the 
normalized,  numerical  results  by  the  exact  values  of  K(S>  for  the  buried 
crack. 

Figures  3 and  4 are  normalized  values  of  K(E)  foi  surface  and 
corner  elliptical  cracks  under  uniform  tension  transverse  to  the  crack  plane. 
Three  elements  combine  to  obtain  K at  any  location:  1)  the  front  face  or 

free  surface  effect  on  K at  the  surface;  2)  the  free  surface  effect  on  K at 
the  deepest  location;  and  3)  the  effect  of  crack  shape  on  K(S).  The  third 
effect  is  partially  contained  in  the  normalization  term;  its  influence  on 
the  normalized  data  in  Figures  3,  4 is  secondary.  In  general,  the  results 
show  that  the  free  surface  magnification  is  around  30-35  percent  at  the  sur- 
face for  both  problems;  the  magnification  at  the  deep  location  is  about  5 
percent  for  surface  cracks  and  12  percent  for  corner  cracks.  Further  numeri- 
cal data  are  given  in  Reference  3. 

Current  numerical  modeling  is  focused  on  evaluating  the  effects  of 
two  other  geometrical  variables:  finite  thickness  and  the  angles  at  the  crack- 

free  surface  intersection.  Correlation  of  the  numerical  results  for  surface 
and  corner  cracks  is  being  generated  by  comparison  of  predicted  fatigue  crack 
growth  to  experimental  results.  Results  to  date  show  good  correlation  at 
both  near-surface  and  deep  crack  locations. 


Kellipse/Kcircle  (NUMERICAL) 
CT:  Kellipse /KcircleCEX ACT) 


FIGURE  2.  VARIATION  OF  STRESS  INTENSITY  FOR  BURIED  CRACK 
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DISCUSSION 

Question;  You  said  ligament  stresses,  if  anything,  go  down  with 
increasing  crack  depth.  How  is  that  so? 

Cruse:  It  is  due  to  the  way  the  crack  sheds  load. 

BIE  computed  stress  intensity  results  are  believed  to  be  too  high. 
I believe  the  experimental  results. 

The  stress  intensity,  by  definition,  does  not  go  to  zero  at  the 
free  surface.  However,  whatever  goes  on  at  the  free  surface  does  not  affect 
the  crack  at  the  root. 

Question:  How  was  K calculated? 


Cruse : From  C.O.D.  next  to  the  crack  tip,  It  was  not  extrapolated 

but  a normalization  procedure  was  used. 

As  the  crack  gets  very  deep,  K seems  to  tail  off  as  C.  W.  Smith 
and  J.  C.  Bell  have  experienced.  However,  I was  not  willing  to  go  beyond 
a/t  of  0.7  because  of  the  dimpling. 

Bell:  With  regard  to  the  back  surface  effect:,  there  are  no  singu- 

larities  at  the  back  surface.  So  what  happens  is  that  as  complications 
occur  on  the  back  surface,  they  influence  the  crack  stress  which  in  turn 
influences  the  front  surface  even  more.  Thus,  in  the  end,  the  front  face 
is  still  the  most  critical  area  for  my  analysis. 

Cruse:  My  analysis  shows  that  the  back  face  stress  intensity 

magnification  factors  are  most  dominatly  tied  to  the  K at  the  front  face. 
Neighboring  effects  are  most  critical  far  elliptical  cracks.  The  most  crude 
local  information  is  obtained  in  the  back  face  model  so  that  I have  trouble 
correlating  deep  crack  K' s wita  experiment. 
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SUBSURFACE  ELLIPTICAL  FLAWS 
by 

A.  S.  Kobayashi,  N.  Polvanich,  A.  F.  Emery, 
and  W.  J.  Love 


The  well-documented  alternating  technique  in  three-dimensional 
fracture  mechanics  was  used  to  obtain  our  recent  results  on  stress  intensity 
factors  of: 

(1)  An  elliptical  crack  near  or  partially  penetrated  through  a 
re-entry  corner  of  a three-quarter  infinite  solid  subjected 
to  a polynomial  distribution  of  stresses; 

(2)  An  elliptical  crack  near  a square  corner  of  a quarter-infinite 
solid  subjected  to  a polynomial  distribution  of  stresses; 

(3)  A large  elliptical  crack  in  a thin  plate  subjected  to  a 
polynomial  distribution  of  stresses. 


Two  iterations  of  the  alternating  technique  yielded  in  ninety 
percent  of  the  crack  surface  residual  surface  tractions  less  than  one  per- 
cent of  the  maximum  prescribed  pressure  for  the  largest  proximity  ratio  of 
b/h  = 0.9.  Residual  surface  tractions  on  the  two  bounding  flat  surfaces  of 
the  solid  were  less  than  one  percent  of  the  maximum  prescribed  pressure. 
Numerical  convergence  for  the  partially  penetrated  elliptical  crack  at  a 
re-entry  corner  was  enhanced  by  prescribing  appropriate  fictitious  pressure 
on  the  penetrated  portion  of  an  elliptical  crack.  The  results  are  presented 
in  terms  of  a stress  intensity  magnification  factor,  which  is  the  actual 
stress  intensity  factor  divided  by  the  corresponding  stress  intensity  factor 
of  an  embedded  illiptical  crack  in  an  infinite  solid,  for  constant  through 
cubic  terms  of  polynomial  crack  pressure  distribution. 

For  an  embedded  elliptical  crack  of  crack  aspect  ratio  of  b/a  = 0.02 
and  a/h  =0.8  and  at  a square  corner,  the  second  bounding  surface  effectively 
increases  the  stress  intensity  magnification  factor,  by  only  0.65  percent 
over  the  corresponding  solution  of  an  embedded  elliptical  crack  in  a semi- 


infinite solid 


(1) 


For  b/a  =0.2  and  a/h  = 0.9,  distinct  influence  of  the 


■a  »<»*•*. 
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second  bounding  surface  is  seen  by  a rapid  increase  in  the  stress  intensity 
magnification  factors  up  to  16  percent.  For  a 20  percent  partially  penetrated 
crack  at  a re-entry  corner,  the  stress  intensity  magnification  factor  in- 


creases to  about  1.11  near  the  free  surface  in  contrast  to  the  1.23  values 

(2) 


observed  in  semi-circular  surface  problems. 


Published  results  on  inner  cracks  in  thermally  shocked  cylinder 


(3) 


unpressurized  inner ^ and  outer  cracks ^ in  pressurized  cylinders,  pres- 


.(4) 


surized  inner  semi-elliptical  cracks  in  pressurized  cylinders as  well  as 
unpublished  results  involving  a comparison  of  calculated  and  experimental 
crack  mouth  openings  in  surface  flawed  plates  subjected  to  tension  and  bending^ 
are  also  presented. 
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DISCUSSION 


F.  W.  Smith:  For  an  ellipse  in  the  middle  of  a plate,  we  also 

showed  that  interaction  is  important.  I agree  with  your  statement  on 
bending  moments  far  away. 

Kobayashi:  I should  have  run  a larger  area  and  not  tried  to 

save  on  computing  time. 

In  the  surface  flaw,  cross  coupling  plays  a large  effect.  The 
model  should  have  been  taken  all  the  way  to  the  grips. 
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STRESS  INTENSITY  FACTORS  FOR  A PRESSURIZED  THICK-WALL  CYLINDER 
WITH  A PART-THROUGH  CIRCULAR  SURFACE  FLAW  - 
COMPLIANCE  CALIBRATION  AND  COLLOCATION  METHOD 


by 

M.  A.  Hussain 


Maggs  Research  Center 
Watervliet  Arsenal 
Watervliet,  New  York 


The  K calibration  at  the  deepest  point  of  a circular  surface  flaw 
in  an  internally  pressurized  thick-walled  cylinder  has  been  determined  by 
a compliance  test  as  well  as  three-dimensional  collocation  method.  Compli- 
ance is  defined  as  the  change  in  internal  volume  of  a cylinder  divided  by 
applied  hydrostatic  pressure  instead  of  usual  elongation/load  definition.  It 
is  shown,  using  linear  theory  of  elasticity  that  the  derivative  Jith  respect 
to  crack  depth  of  internal  and  external  volume  changes  are  identical.  This 
permits  the  use  of  external  strain  measurements  to  compute  the  stress  intensity 
factors.  Periodic  cubic  splines  are  used  to  interpolate/approximate  the 
strain  data  as  well  as  to  compute  external  volume  change  as  a function  of 
, crack  depth. 

With  a starter  surface  notch  of  semi-circular  cross-section,  an 
extensive  set  of  tests  were  performed  on  a cylinder  with  7.1  inch  smooth  bore 
and  14.25  inch  outside  diameter.  A fatigue  crack  was  grown  from  this  starter 
notch  until  it  reached  each  successive  desired  depth.  The  circumferential 
strains  were  then  read  on  outside  circumference  at  14  angular  locations.  The 
stress  intensity  factors  were  computed  by  compliance  technique. 

The  results  at  the  deepest  points  compared  well  with  exterior 

collocation  method.  In  this  method  we  represent  the  shape  of  such  a crack 

by  a segment  of  a circle.  The  objective  is  to  construct  a set  of  functions 

such  that  the  most  stringent  boundary  conditions  in  the  plane  of  crack  are 

satisfied  exactly.  These  conditions  are:  the  normal  displacement  u^  is  zero 

outside  the  circle  c-'d  the  normal  stress  o is  zero  inside  the  circle.  For 

z 

a symmetric  problem,  the  shearing  stresses  vanish  in  the  entire  crack  plane. 
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Using  integral  forms  of  harmonic  functions,  the  above  conditions 
are  satisfied  and  the  rest  of  the  boundary  conditions  exterior  to  the  plane 
of  the  crack  are  then  satisfied  by  collocation  in  the  least  square  sense. 

This  is  analogous  to  the  two-dimensional  problem  where  William's 
stress  functions  have  been  used  successfully.  However,  as  opposed  to  two- 
dimensions,  the  three-dimensional  surface  flaw  problems  require  an  additional 
condition  to  be  satisfied;  namely  the  boundary  traction  must  vanish  at  points 
of  crack-free  surface  intersection.  In  our  analysis  this  requirement  yields 
a zero  stress  intensity  factor  at  such  points.  However,  the  singularity  at 
crack-free  surface  intersection  points  has  yet  to  be  investigated. 


DISCUSSION 


Question:  Do  the  results  represent  one  point  or  an  average  of 

the  K's  taken  around  the  crack  front? 

Husain:  The  results  are  taken  only  at  z = 0. 

Cruse:  Are  not  the  strain  gages  on  the  mid  crack  line  measuring 

an  average  SIF  while  your  analytical  data  is  single  valued? 

Husain:  The  strain  gages  are  not  measuring  an  average. 

F,  W.  Smith:  Olavi  in  about  1966  published  results  lor  a circular 
crack  problem  for  which  the  surface  intersects  off-center  and  obtained 
different  results. 
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ON  THE  THREE-DIMENSIONAL  THEORY  OF  FRACTURE 


by 

E.  S.  Folias 

University  of  Utah 
Salt  Lake  City,  Utah 


The  speaker  discussed  his  recent  theoretical  results  on  the  linear 
elastic,  three-dimensional,  stress  distribution  in  a plate  of  finite  thick- 
ness and  containing  a plane  (rectangular),  through-the-thickness , crack  and 
under  the  action  of  Mode  I loading. 

Emphasis  was  placed  on  the  physical  interpretation  of  the  results 
and  their  implications  to  fracture. 

Finally,  the  results  were  compared  with  currently  existing  in 
the  literature  theoretical  and  experimental  evidence  and  attention  was  drawn 
to  some  possible  pitfalls. 
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DISCUSSION 


Cruse:  It  is  my  contention  that  all  analytical  studies  break 

down  at  the  surface  corner  of  the  through  crack  because  the  discontinuity 
of  the  behavior  is  built  in.  Thus,  you  cannot  get  unique  results. 

Folias:  Sternburg  & Sadowsky  get  K going  down  at  the  corner  for 

different  geometries  than  I have,  but  these  are  not  comparable. 

Sih:  How  do  you  define  uniqueness? 

Folias:  Wilcox's  proof  is  short  and  I can  send  it  to  anyone 

interested.  He  has  worked  extensively  with  the  defraction  equation  and  was 
the  first  one  to  prove  uniqueness  for  this  problem. 
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SOME  UNSOLVED  SINGULARITY  PROBLEMS 
by 

M.  L.  Williams 

University  of  Pittsburg 
Pittsburg,  Pennsylvania 


A key  to  understanding  fracture  problems  is  to  understand  the 
nature  of  the  singularity  existing  in  the  field.  Consideration  of  some  of 
the  following  unsolved  singularity  problems  could  lead  to  increased  under- 
standing of  the  singularities  in  3-d  fracture. 

(1)  A conical  hole  in  a half  space  into  which  a rigid  sphere 
is  pressed.  The  problem  of  whether  there  is  a singularity 
at  the  tip  of  the  conical  hole  has  not  been  resolved  math- 
ematically although  there  appears  to  be  no  singularity. 

(2)  The  mathematically  important  question  of  completeness  for 
the  stress  functions  at  a V-notch  in  a plane  stress  speci- 
men has  not  been  solved.  This  question  also  has  not  been 
answered  for  fracture  in  adhesive  bonded  joints. 

(3)  The  problem  of  the  oscillation  of  the  singularity  near 
the  crack  tip  for  mixed  media  (adhesive)  problems  makes 
the  choice  of  the  proper  size  of  the  finite-elements  in 
this  crack  problem  very  difficult. 

(4)  A problem  posed  by  Rongved  of  a square  or  rectangular  block 
bonded  to  a naif  space  with  a normal  tensile  load  applied 
to  the  block.  What  is  the  singularity  at  the  corners  where 
both  the  geometrical  and  mathematical  boundary  conditions 
are  not  smooth? 

(5)  In  a center  cracked  plate  of  finite  width  what  is  he  behavior 
as  the  crack  nears  the  edge  (The  last  fiber  problem)? 

(6)  What  is  the  behavior  in  a clamped  free  right  angled  corner 

as  the  geometry  passes  L'niough  the  transition  from  no  initial 
crack  to  a finite  length  crack  along  the  clamped  edge. 
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FUNDAMENTAL  STUDY  OF  CRACK  INITIATION  AND  PROPAGATION 

by 

Mark  L.  Wilkins 

Lawrence  Livermore  Laboratory 
Livermore,  California 


The  program  was  started  at  the  Lawrence  Livermore  Laboratory  to 
predict  the  capabilities  and  potential  failure  of  a nuclear  reactor  pres- 
sure vessel.  A coupled  calculational  and  experimental  approach  is  used  to 
predict  initial  failure.  The  work  includes  calibrating  a fracture  model 
in  two-  and  three-dimensional  time  dependent  computer  programs  used  to 
simulate  engineering  tests  of  material.  Examples  show  application  of  the 
three-dimension  HEMP  program  to  calculate  the  critical  state  of  failure  of 
a flat  tensile  specimen  and  the  calculated  stress  field  around  a through - 
the-thickness  crack  in  a pipe. 
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APPLICATION  OF  AN  INFLUENCE  FUNCTION  METHOD  FOR 
THREE-DIMENSIONAL  ELASTIC  ANALYSIS  OF  CRACKS 

by 

Philip  M.  Besuner 

Failure  Analysis  Associates 
Palo  Alto,  California 


A weight  or  influence  function  (IF)  method  has  been  extended  to 
allow  accurate  and  extremely  inexpensive  calculations  of  elastic  stress 
intensity  factors  in  three-dimensional  crack  problems.  The  stress  intensity 
factors  are  not  local  values,  K(s)  at  some  position  on  the  crack  front  s, 
but  rather  are  local  root-mean-square  values,  K.  The  K calculations  are  in- 
expensive numerical  integrations  of  the  product  of  predetermined  influence 
functions  and  uncracked  stress  functions  (i.e.,  the  stress  field  at  the  crack 
locus  in  the  uncracked  solid).  The  method  allows  any  number,  n,  of  these 
rms  stress  intensity  factor  values,  K^;  i = 1,  n to  be  calculated  around  the 
crack  front.  In  the  limit  of  asymptotically  large  n,  it  can  be  shown  that  K(s) 
approaches  the  local  value  K(s).  In  this  method,  the  crack  is  described  by 
any  number  of  scalar  dimensions  a^;  i = 1,  n or  degrees  of  freedom  (DOF), 
and  each  corresponds  to  the  strain  energy  release  rate  obtained  by 
perturbation  of  a single  DOF  while  holding  all  other  DOFs  constant. 

The  advantages  of  the  methou  derive  from  the  fact  that  influence 
functions  may  be  calculated  through  the  crack  opening  displacements  deter- 
mined analytically  or  experimentally  for  the  desired  crack  geometry  and  dis- 
placement constraints  under  very  simple  loading,  such  as  uniform  pressure 
inside  the  crack.  These  influence  functions  may  then  be  used  to  compute  K^ 
for  a general  bivariate  uncracked  stress  field,  such  as  o^(x,y).  Thus,  the 
method  has  the  important  advantage  that  neither  error  nor  cost  of  solution 
is  affected  by  the  complexity  of  the  known  uncracked  stress  field. 

The  author  presented  justification  of  the  use  of  K and  summarized 
many  successful  engineering  applications  of  the  IF  method  in  which  one  to 
four  DOFs  have  been  used  to  model  the  three-dimensional  crack.  Most  of  the 
applications  have  involved  elliptical  cracks  with  multiple  degrees  of  freedom 
under  fatigue  loading.  By  computing  the  growth  rate  of  each  dimension  a_^ 


from  its  associated  K^,  the  method  has  the  capability  to  predict  non-self- 
similar fatigue  crack  growth,  as  is  indicated  in  Figure  1,  which  illustrates 
a two  DOF  quarter-elliptical  corner  crack. 

An  application  to  the  analysis  of  boiling  water  reactor  feedwater 
nozzle  corner  cracks  was  outlined.  Both  the  IF  methoc.  and  a readily  avail- 
able, but  not  optimum,  finite-element  (FE)  method  were  applied  to  several 
three-dimensional  crack  analysis  test  cases.  The  test  cases  had  been  chosen 
for  their  similarity  to  a nozzle  crack  analysis  being  funded  by  the  Electric 
Power  Research  Institute  and  for  the  availability  of  an  accurate  published 
result  obtained  from  some  recognized  third  method  of  solution.  Results  were 
given  which  summarize  bcMi  accuracy  and  computer  cost  of  the  two  methods  for 
each  of  the  selected  cases.  The  IF  method  was  demonstrated  to  be  superior 
from  both  an  accuracy  and  cost  viewpoint.  The  FE  method,  as  applied,  was 
more  than  1,000  times  more  costly  than  the  IF  method  and  had  average  errors 
of  3-5  percent  *>nd  peak  errors  of  8 percent,  as  compared  to  less  than  2 per- 
cent error  for  the  IF  method.  On  the  basis  of  these  results,  the  IF  method 
has  been  chosen  by  Electric  Power  Research  Institute  to  perform  all  stress 
intensity  factor  evaluations  of  nozzle  flaws  in  the  remaining  phases  of  a 
current  industry  study. 


Local  Value  3 K (0,  a^) 


FIGURE  1.  TWO-DEGREE-OF-FREEDOM  ELLIPTICAL  CORNER  CRACK  MODEL 
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ELLIPTICAL  CRACK,  SURFACE  FLAW,  AND  FLAWS  AT  FASTENER  HOLES 

by 

Ramesh  C.  Shah 

Boeing  Aerospace  Company 
Seattle,  Washington 


This  paper  summarizes  the  work  done  on  evaluations  of  elastic 
stress  intensity  factors  of  three-dimensional  crack  problems  such  as  embedded 
elliptical  cracks,  semi-elliptical  surface  flaws  and  part-through  cracks 
originating  at  fastener  holes. 

Stress  intensity  factors  for  an  embedded  elliptical  crack  approach- 
ing the  free  surface  of  a semi- inf inite  solid  and  subjected  to  a uniform 
pressure,  and/or  a linearly  varying  pressure  are  determined  in  a nondimen- 

sional  form  as  a function  of  position  around  the  crack  periphery,  crack 

(l  2) 

aspect  ratio,  and  crack  distance  from  the  free  surface  ’ . Stress  inten- 

sity factors  are  determined  numerically  using  an  alternating  technique  with 

(3) 

two  solutions  . The  first  solution  involves  an  elliptical  crack  in  a solid 
and  subjected  to  normal  loading  expressible  in  a polynomial  of  x and  y^\ 

The  second  solution  involves  stresses  in  the  half  spac;  due  to  prescribed 
normal  and  shear  stresses  on  the  surface^.  Effect  of  the  Poisson's  ratio 
of  the  material  on  these  stress  intensity  factors  is  also  investigated  . 
Nondimensional  stress  intensity  factors  for  the  elliptical  crack  in  a semi- 
infinite solid  subjected  to  uniform  tension  are  compared  with  those  derived 
by  Nisitani  and  Murakami ^ and  agreement  between  the  two  is  excellent.  Non- 
dimensional  stress  intensity  factors  for  an  elliptical  crack  in  a finite 
thickness  plate  subjected  to  tensile  and  bending  loadings  are  determined  from 
the  above  solutions  by  assuming  negligible  coupling  effects  between  the  two 

/ O \ 

free  surfaces v . These  results  compare  very  well  with  those  evaluated  by 

(9) 

Smith  by  the  alternating  technique  with  both  stress  free  surfaces  present. 

Nondimensional  stress  intensity  factors  at  the  maximum  crack  depth 

of  a semi-elliptical  surface  crack  in  a finite  thickness  plate  subjected  to 

tensile  and  bending  loadings  are  estimated  for  various  crack  aspect  ratios 

0 2) 

and  crack  depth  to  plate  thickness  ratios  ’ . These  nondimensional  stress 


WW. 
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intensity  factors  are  compared  with  those  obtained  by  other  investigators 


(7,10-15) 


.(17) 


and  some  of  these  results ^»-'-0,15,16)  are  usg£j  tQ  compUte  fracture  toughness 
values  of  static  fracture  tests  of  aluminum  and  titanium  surface  flawed  speci- 
mens in  tensiou^^.  Results  show  that  stress  intensity  solutions  of  Refer- 
ences^’^’^^  produce  reasonably  constant  values  of  fracture  toughness. 

A procedure  is  formulated  to  derive  approximate  stress  intensity 
factors  for  two  embedded  semi-elliptical  or  through-the-thickness  cracks 
originating  at  open  holes  or  at  loaded  or  unloaded  fastener  filled  holes* 

The  procedure  involves  obtaining  the  stress  distribution  at  locations  of  two 
semi-elliptical  cracks  in  an  uncracked  solid  with  the  fastener  hole  for 
appropriate  loading  conditions  and  then  pressurizing  an  elliptical  crack  in 
a solid  without  the  fastener  hole  with  this  stress  distribution.  Nondimen- 
sional  stress  intensity  factors  are  derived  for  through  cracks  at  loaded  close 
tolerance  fastener  filled  holes and  at  interference  fit  fastener  filled 
holes'  , for  semi-elliptical  cracks  at  open  holes  and  at  loaded  and  unloaded 

close  tolerance  fastener  filled  holes^\  and  for  semi-circular  cracks  at 

(19) 

interference  fit  fastener  filled  holes  . When  the  problem  of  two  semi- 
elliptical cracks  at  an  open  hole  is  reduced  to  a two-dimensional  crack  problem 

of  t .-ough-the-thickness  cracks  by  letting  a /c->=%  stress  intensity  factors 

(20) 

agree  very  well  with  those  given  by  Bowie  . Solutions  of  two  embedded  semi- 
elliptical cracks  at  fastener  holes  in  a solid  are  made  suitable  to  quarter- 
elliptical  corner  cracks  at  fastener  holes  in  a finite  thickness  plate  by 
applying  available  appropriate  free  surface  correction  factors.  These  results 

were  applied  to  static  fracture  and  fatigue  crack  propagation  tests  containing 

(18 

corner  cracks  at  fastener  holes  for  various  fastener  and  loading  conditions'  ’ 
18  21  22) 

’ ’ , and  agreement  between  computed  and  actual  results  is  very  good. 
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AN  ASSESSMENT  OF  NEAR  TIP  MODELING  FOR 
2D  AND  3D  CRACK  PROBLEMS 

by 

C.  F.  Shih 

General  Electric  Company 
Schenectady,  New  York 

There  are  several  approaches  to  modeling  the  singular  field  at 
the  crack  tip  via  the  displacement  approach.  For  purposes  of  discussion, 
these  may  be  grouped  into  three  categories. 

In  the  first  category,  the  "entire"  singularity,  i.e.,  the  radial 
and  circumferential  dependence  of  the  singular  fields  are  embedded  through 
the  use  of  special  elements  or  singular  elements.  These  elements  were  used 
by  Wilson,  Hilton,  Hutchinson,  Shih,  etc.  More  recently,  Benzley  introduced 
the  "enriched"  element,  which  is  an  isoparametric  element  that  contains  both 
the  usual  interpolation  functions  as  well  as  the  singular  functions. 

The  next  category  of  special  elements  enforces  the  radial  depen- 
dence ^f  the  singular  field  but  allows  tlv’  circumferential  dependence  to  be 
determined  by  the  finite-element  method.  Such  elements  were  used  by  Tracy, 
Levy,  Rice,  etc. 

In  the  last  category,  is  the  "quarter  point"  side  node  isopara- 
metric element  which  gives  the  square  root  singularity  along  the  edge  of  the 
"misplaced"  side-nodes.  These  elements  have  been  employed  by  Henshell  and 
Shaw,  Barsoum,  Bloom,  etc. 

The  above  summarizes  some  of  the  more  well-known  elements  that 
are  employed  to  model  near  tip  fields.  Each  category  of  elements  have  their 
advantages  and  disadvantages  and  these  were  discussed.  The  accuracy  of  the 
above  elements  will  be  compared  for  several  standard  crack  configurations. 
Some  of  the  author's  preliminary  work  on  modeling  crack  growth  for  through 
cracks  were  presented. 
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HYBRID  MODELS  FOR  THREE-DIMENSIONAL 
CRACK  ELEMENTS 

by 

Theodore  H.  H.  Plan 

Massachusetts  Institute  of  Technology 
Cambridge,  Massachusetts 


This  presentation  will  first  review  the  various  modified  varia- 
tional principles  in  solid  mechanics  and  the  corresponding  finite-element 
hybrid  model  for  the  development  of  special  singular  elements  for  linear 
fracture  mechanics.  The  basic  schemes  for  three  different  hybrid  models 
are : 


(1)  A scheme  which  is  based  on  the  assumed  equilibrating  stress 
field  which  also  satisfies  the  compatibility  condition  inside 
the  element,  and  on  independently  assumed  boundary  displace- 
ments. Such  model  may  be  interpreted  as  either  a hybrid 
stress  model  or  a hybrid  displacement  model. 

(2)  A scheme  which  is  based  on  assumed  equilibrating  stress 
field  inside  the  element  and  independently  assumed  boundary 
displacements.  Such  model  is  a hybrid  stress  model. 

(3)  A scheme  which  is  based  on  assumed  displacement  field  inside 
the  element  and  independently  assumed  boundary  displacements. 
Such  a model  is  a hybrid  displacement  model. 

The  applications  of  these  three  schemes  to  two-dimensional  crack 
elements  have  been  reported  by  Tong,  Pian  and  Lasry  (1973);  by  Pian,  Tong 
and  Luk  (1971);  and  by  Atluri,  Kobayashi  and  Nakagaki  (1974),  respectively. 

It  appears  that  the  most  desirable  element  is  one  formulated  by  the  first 
scheme  for  which  an  imbedded  crack  is  included  in  the  element  and  only  the 
nodal  displacements  are  left  as  unknowns.  In  the  implementation  of  such  an 
element,  only  boundary  integrations  are  called  for.  Since  the  boundary  does 
nut  xueluue  tut:  clack  lip,  the  integration  does  not  involve  any  singular  terme 
The  application  of  the  other  two  schemes,  however,  will  result  in  the  use  of 
several  elements  around  the  tip  of  the  crack.  \t  is  then  needed  to  introduce 
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the  stress  intensity  factors  as  unknowns  in  addition  to  the  nodal  displace- 
ments. Furthermore,  the  evaluation  of  the  stiffness  matrices  in  these  cases 
involves  boundary  integration  with  singular  terms  and  special  treatments  are 
needed  in  the  numerical  integration  process. 

The  main  discussion  will  be  on  the  extension  of  these  three  models 
to  three-dimensional  crack  elements.  Here,  however,  the  first  two  schemes 
may  become  difficult  because  of  their  requirement  of  assumed  equilibrating 
stresses.  It  is  particularly  so  when  the  problem  involves  a curved  line  as 
its  crack  front.  However,  for  an  element  with  a straight  crack  front,  it  is 
not  difficult  to  construct  equilibrating  stresses  that  will  provide  the 
correct  near  field  singular  behavior. 

In  all  the  formulations  of  hybrid  elements  for  three-dimensional 
crack  problems,  it  is  generally  necessary  to  evaluate  some  surface  integrals 
that  involve  singular  terms.  Also  because  of  the  difficulty  in  constructing 
the  singular  solutions  to  satisfy  both  the  equilibrium  and  compatibility 
conditions,  it  is  further  complicated  by  the  need  for  performing  volume  inte- 
grations that  involve  singular  terms. 

The  presentation  included  a brief  report  of  results  of  crack  tip 
stress  distributions  of  graphite  fiber  reinforced  composite  obtained  by  using 
hybrid  stress  elements,  which  do  not  include  stress  singularities.  Those 
include  the  use  of  laminated  plate  elements  which  take  into  account  inter- 
laminate shear  stresses  and  three-dimensional  : >lid  elements.  A key  feature 
for  the  success  of  all  these  elements  is  the  inclusion  of  the  stress  free 
conditions  in  the  element  development.  Such  conditions  can  be  easily  accom- 
plished by  the  assumed  stress  hybrid  method.  The  analysis  includes  the 
modeling  of  a 90/0/0/90  laminate  with  subcritical  splits  propagating  parallel 
to  the  fibers  of  each  ply.  For  laminated  composites  the  near  field  stress 
distribution  includes  singular  behavior  due  to  t he  presence  of  the  c k and 
the  effect  of  interlaminar  stresses  at  the  free  edge.  Thus,  it  is  uncertain 
whether  special  singular  elements  can  be  developed  for  such  applications. 


wm 


1 


1! 

I 

I 

I 

i 

I 

I 

I 


f 


) 


*t 

I 

AM 


84 

THE  FINITE-ELEMENT  ALTERNATING  METHOD  FOR 
ANALYSIS  OF  COMPLEX  THREE-DIMENSIONAL  CRACK  PROBLEMS 

by 


F.  W.  Smith 

Colorado  State  University 
Ft.  Collins,  Colorado 


A new  procedure  for  analysis  of  three-dimensional  crack  problems 
is  described  which  makes  the  use  of  the  f inite-e] tment  method  in  conjunction 
with  the  alternating  method.  The  procedure  combines  the  best  features  of 
closed  form  elasticity  solutions  for  crack  problems  and  the  finite-element 
method  to  provide  an  analysis  tool  which  is  flexible,  relatively  inexpensive 
and  appears  to  be  reasonably  accurate. 

The  alternating  procedure  operates  by  applying  stresses  to  the 
crack  surface  and  using  the  known  crack  solutions  to  compute  stresses  in  a 
semi- inf inite  solid  at  the  locations  of  the  surfaces  of  the  solid  in  which 
the  crack  is  assumed  to  reside.  The  opposite  of  these  stresses  is  then 
applied  to  the  surface  of  the  body  in  question,  freeing  the  surface  of  stress, 
and  the  finite-element  method  is  used  to  compute  resulting  stresses  at  the 
location  of  the  crack  surface.  The  crack  surface  stresses  are  again  removed 
as  before,  and  the  process  is  repeated  until  further  contributions  are  negli- 
gible. The  results  of  all  iterations  are  then  superirr  osed  to  obtain  the 
crack  opening  displacements  and  the  stress  intensity  factors. 

Check  cases  have  been  analyzed  which  provide  a comparison  of 
results  from  this  procedure  with  accepted  two-dimensional  results.  Agreement 
in  these  cases  is  within  three  percent. 

A series  of  results  is  presented  for  single  and  double,  corner 
and  embedded  surface  flaws  near  holes  in  flat  plates  subjected  to  uniform 
tension.  These  results  are  quantitatively  compared  with  experimental  results 
based  on  stress  freezing  and  on  fatigue  studies  of  cracks  having  similar 
geometries.  The  results  are  qualitatively  compared  with  high  speed  motion 
pictures  taken  of  the  initial  slow  growth  occurring  in  plexiglass  having 
cracks  near  a hole.  The  strengths  of  the  procedure  presented  are  as  follows: 
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It  is  based  on  continuum  crack  solutions  and  gives  the  correct 
stress  behavior  at  the  crack  tip. 

It  is  flexible  with  respect  to  problem  geometry. 

It  is  cost  effective  in  terms  of  computing  costs.  After  set  up, 
a solution  uses  220  seconds  of  CPU  time. 

It  is  accurate  in  comparison  with  known  solutions. 


The  weaknesses  are: 


It  uses  a polynominal  fitting  scheme  to  deal  with  stresses  on 
the  crack  surface  which  may  cause  errors  in  problems  other  than 
the  check  cases. 

It  does  not  deal  with  the  details  of  the  problem  of  the  crack 
front  intersection  with  a free  surface. 
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SESSION  V 

PRACTICAL  PROBLEMS  OF  INTEREST  TO  GOVERNMENT  AGENCIES 


Session  V consisted  of  presentations  by  representatives  of  five 
Government  agencies.  These  representives  described  various  three-dimen- 
sional fracture  problems  of  high  interest  to  their  respective  agencies. 

The  government  agencies  who  were  represented  in  session  were 

Army  (Watervliet,  AMMRC) 

Air  Force  (AFFDL , AFML,  AFRPL,  AFOSR) 

NASA  (Lewis  Research  Center) 

Department  of  Transportation  (TSC) 

Oak  RLdge  National  Laboratory  (ERDA) 

In  general,  these  discussions  were  quite  brief.  A summary  of  each 
representative’s  presentation  is  given  below. 

Army  (Watervliet,  M.  A.  Hussain) 

Geometries  of  flaws  of  interest  at  Watervliet  are  semi-circular  or 
elliptical  surface  cracks  that  occur  in  gun  tubes  which  can  coalesce  into 
larger  frontal  cracks.  Prediction  of  the  rate  of  formation  of  these  cracks 
is  complicated  by  the  presence  of  residual  stresses  which  are  developed 
during  the  gun  tubes  manufacture  (an  auto  frett  ,'ng  process). 

To  a question  on  the  dependence  of  crack  propagation  in  gun  tubes 
to  strain  rates  the  answer  was  given  that  it  is  not  presently  thought  to  be 
an  important  consideration.  It  was  brought  out  that  erosion  of  the  inner 
gur  tube  surface  is  significant  and  of  concern  in  itself  but  since  its  occur- 
rence is  largest  at  the  biggest  end  of  the  gun  tubes  it  is  not  considered  to 
be  a significant  structural  phenomena. 

Army  (AMMRC,  D.  M.  Tracey) 

AMMRC  is  concerned  about  structural  and  material  modeling  of  cracks 
in  various  hardware  components  such  as  guns,  missiles,  jeeps,  etc.  These 
cracks  typically  result  from  fabrication,  fatigue,  corrosion  or  wear. 
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It  is  AMMRC’s  expressed  op  “lion  that  the  effects  upon  crack  growth 
that  have  not  been  adequately  addressed  in  the  literature  are  those  due  to 
inertial  loading  (such  as  in  impact  problems)  and  residual  stress.  Further- 
more, because  of  the  three-dimensionality  of  the  state  of  stress  in  many 
problems,  mixed  mode  crack  propagation  can  not  be  neglected.  In  view  of 
greater  use  of  composite  materials,  the  effects  of  anisotropy  and  nonhomo- 
geneities, (layers,  coatings)  are  very  real  problems  and  in  need  of  treatment. 
It  was  expressed  that  some  of  our  efforts  would  be  well  spent  if  we  would 
attempt  solutions  to  these  new  aspects  of  crack  problems  rather  than  make 
small  improvements  on  present  ones. 

It  was  also  expressed  that  there  exists  a high  degree  of  personali- 
zation in  our  present  day  publications  and  that  we  have  not  made  specialized 
tools  (computer  codes)  generally  available.  These  tools  are  often  undocu- 
mented. We  should  be  more  conscious  of  this  lack  of  documentation  of  soft- 
ware and  seek  support  to  document  each  of  our  individual  capabilities  in 
order  that  others  can  make  full  advantage  of  them. 


Air  Force  (AFFDL,  H.  A.  Wood) 

The  new  design  phiJosophy  instituted  by  the  Air  Force  involves 
damage  tolerance  analysis,  i.e.,  it  is  assumed  that  the  structure  contains 
flaws  in  its  early  life.  It  is  felt  that  to  insure  safe  operation  of  air- 
craft is  to  insure  that  propagation  of  these  flaws  does  not  become  catas- 
trophic in  some  time  period.  This  time  period  could  be  the  optimum  service 
life  or  the  service  interval  between  inspections.  Specification  of  new 
designs  re.quire  that  the  engineer  design  the  structure  for  residual  strength 
and  safe  life.  The  procedure  is  to  assume  that  at  time  zero  there  are  flaws 
in  the  structure  and  that  the  size  of  the  flaw  is  dependent  upon  one’s  ability 
to  inspect  the  structure.  Two  requirements  for  failure  are  that 

(1)  Crack  growth  life  exceeds  some  design  value  or 

(2)  Residual  strength  equals  or  falls  short  of  some  minimum 
value. 

During  the  damage  tolerance  analysis,  it  is  assumed  that  the  flaws 
exist  in  the  highest  stress  field  areas  and  that  they  are  orientated  in  the 
worst  possible  direction. 
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A most  prevalent  flaw  geometry  -“s  one  that  manifests  itself  at 
holes  in  wing  spars. 

Areas  of  general  interest  to  the  Air  Force  are  the  effects  upon 
cra-'k  growth  rate  of  material  variations,  load  history  effects  (such  as 
variable  amplitude  and  real  time  effects),  and  crack  tip  placticity.  The 
general  three-dimensional  geometries  of  interest  are  corner  cracks,  thru- 
and  part-thru  cracks,  cracks  at  stress  concentrations  and  the  effects  of 
the  front  and  back  surface  upon  crack  propagation.  Study  of  these  flaws  is 
made  more  difficult  by  the  complex  loadings  that  can  occur  (such  as  during 
load  transfer  in  bolted  joints  and  interference  fit  fasteners). 

The  basis  of  present  criteria  for  crack  growth  are  totally  elastic 
and  are,  therefore,  subject  to  question.  In  problems  involving  three-dimen- 
sional flaws,  such  as  embedded  cracks,  correlation  of  the  variation  of  K 
around  the  flaw  to  crack  growth  data  is  of  great  interest.  It  is  felt  that 
there  is  a need  for  a broad  new  basis  for  experimental  techniques  to  verify 
K solutions  for  these  problems. 

To  the  question  of  whether  the  Air  Force  assumes  that  flaws  exist 
at  all  holes,  the  answer  was  given  that  only  one  primary  damage  site  is 
assumed.  For  example,  on  a 10-foot  wing  spar  only  one  primary  flaw  is 
assumed.  If  necessary,  depending  upon  the  flaw  size,  the  structure  may  be 
analyzed  assuming  that  the  flaw  can  influence  the  local  load  transfer  in 
the  structure. 

Air  Force  (AFML,  T.  Nicholas) 

The  Air  Force  Materials  Laboratory  has  also  adopted  the  damage 
tolerant  design  philosophy  for  the  design  of  engine  components.  A major 
difference  in  the  nroblems  encountered  in  engine  components  that  are  not 
encountered  in  air  frame  components  is  in  the  type  of  loading.  In  engine 
components,  both  mechanical  and  thermal  cycles  are  experienced,  during  which 
time,  time-dependent  plastic  deformations  and  creep  can  create  strains  as 
large  as  cne  to  two  percent.  Under  this  load  environment,  cracks  are 
believed  to  propagate  under  iow  cycle  fatigue. 

Thus  far,  the  approach  at  AFML  has  been  to  use  a linear  fracture 
mechanics  approach  to  these  complex  problems,  i.e.,  it  is  assumed  tha:  AK 
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and  the  R ratio,  0m£n/°max>  are  important  parameters  (temperature  is 
also  considered).  The  validity  of  this  approach,  however,  is  very  much 
in  question. 

The  geometry  of  problems  of  interest  to  AFML  (like  those  in  the 
air  frame  structures  of  AFFDL)  are  cracks  that  originate  at  holes.  Of  great 
concern  is  what  is  the  critical  crack  length  and  how  long  does  it  take  for 
a crack  to  reach  critical  size? 

The  analysis  technique  in  use  at  AFML  as  at  AFFDL  is  linear  frac- 
ture mechanics.  However,  it  is  AFML's  opinion  that  a different  approach  is 
necessary  for  their  problems. 

In  the  discussions,  it  was  pointed  out  that  tor  the  design  of 
engine  components  there  is  a need  for  greater  emphasis  on  developing  the 
criteria  for  crack  growth  and  modeling  of  crack  growth.  Present  data  in 
terms  of  AK  and  da/dN  is  not  sufficient. 

Air  Force  (AFRPL , R.  Peeters) 

Staff  at  aFRPL  are  presently  funding  studies  in  failure  mechanisms, 
including  initiation,  velocity  of  propagation,  and  stability  of  cracks  in 
solid  propellant  rocket  motors.  These  problems  are  complicated  by  the  fart 
that  propellants  often  exhibit  viscoelastic  behavior.  The  analysis  of  these 
cracks  is  further  complicated  by  the  loading  environment. 

Ore  example  of  the  complex  nature  of  these  problems  is  that  of 
cracks  under  igritior  conditions  (i.e.,  burning  cracks).  In  this  problem, 
the  crack  experiences  very  large  thermal  and  pressure  loads  on  its  faces. 

Other  problem  areas  briefly  mentioned  were  related  to 

• Development  of  a three-dimensional  finite-element  computer 
code  with  specialized  capability  for  the  analysis  of  cracks. 

This  code  will  contain  a displacement  hybrid  finite-element 
capability  for  both  homogeneous  and  bi-material  problems. 

• Development  of  a constitutive  theory  for  materials  with 
process  induced  bonding  states  (i.e.,  nonlinear  viscoelastic). 

• The  study  of  failures  at  the  case/liner  bond  region  of  solid 
propellant  rocket  motors. 
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Air  Force  (AFOSR,  W.  J.  Walker) 


AFOSR  is  striving  to  develop  the  capability  of  life  prediction 
througn  damage  tolerance  analyses.  Current  efforts  include  experimental 
and  analytical  mechanical  studies  and  metallurgical  studies  to  determine 
the  proper  measure  of  crack  growth  characterization,  i.e.,  AK,  energy 
release  rate,  etc. 

In  these  studies,  it  is  an  accepted  fact  that  errors  will  be 
present,  therefore,  it  rs  believed  that  extremely  exact  calculations  are 
not  necessary.  What  is  necessary,  however,  is  that  the  basic  phenomena 
are  modelled  correctly.  One  such  phenomenon  of  interest  to  AFOSR  is  free 
edge  behavior  and  in  determining  if  free  edge  behavior  is  a driving  force 
in  crack  growth. 

Other  problems  that  were  said  to  be  of  interest  to  AFOSR  per- 
tained to 

• The  back  face  problem  and  the  effect  of  methods  of  its 
analysis  upon  design  life  production 

• The  effect  of  plasticity  upon  crack  growth 

• Nonhomogeneous  material  behavior  as  encountered  in 
bonded  layer  construction 

• Biaxial  loading  effects 

• The  presence  of  multiple  cracks  as  in  ceramics. 
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NASA  (Lewis  Research  Center,  C.  C.  Chamis) 


A great  concern  was  expressed  for  the  need  of  sound  structural 
analysis  methods  to  determine  the  load  distribution  and  transfer  through  the 
structure  before  an  analysis  of  the  crack  is  attempted. 

An  interest  was  expressed  in 

• Embedded  and  thru  crack  geometries 

• Mixed  mode  behavior  (both  for  analytical  methods  and 
experimental  behavior) 

• Effects  of  residual  and  cyclic  stresses. 

It  was  felt  that  benchmark  problems  worthy  of  consideration  are 
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0 A surface  crack  in  a finite  thickness  and  width  plate  under 
tension  and  shear 

a A surface  crack  in  a thick  wall  cylinder  or  tube. 

Department  of  Transportation  (TSC , D.  McConnell,  E.  Savage) 

TSC,  Transportation  Systems  Center,  is  a part  of  the  Department 
of  Transportation.  The  work  at  TSC  is  supported  by  a number  of  federal 
agencies  including  FRA  and  FAA. 

For  FRA,  Federal  Railway  Administration,  TSC  is  involved  in  several 
problem  areas  involving  failure  of  components  necessary  to  rail  transportation. 
These  incl ude 

• Wheel  failures  on  rail  cars.  This  is  a thermally  driven 
problem  caused  by  braking. 

• Tank  car  failures.  These  are  in  the  form  of  puncture  or 
penetration  caused  by  couplers  during  collisions  cf  derail- 
ments. Another  form  of  failure  of  interest  is  related  to 
fatigue  crack  propagation  of  flaws  in  weldments  at  the  head/ 
body  intersection. 

• Rail  failures.  This  is  the  most  active  area  at  present  and 
a more  detailed  description  follows. 

TSC's  objectives  are  to  understand  the  areas  of  the  rail  subject 
to  rail  failures,  determine  the  population  of  failures,  determine  the  rate 
of  crack  growth,  and  the  requirements  for  early  detection.  The  approach 
th'  has  been  taken  is  to  attempt  application  of  present  day  technology  to 
understanding  the  behavior  of  cracks  in  rails  rather  than  to  develop  new 
techniques  for  their  analysis. 

Due  to  steadily  increasing  wheel  loads,  as  a result  of  increase 
in  car  weight,  and  decrease  in  track  maintenance,  accident  damage  trends 
have  steadily  increased  during  recent  years  (Figure  1).  This  damage  is 
primarily  attributed  to  rail  failures,  (Figure  2). 

Rail  failures  are  primarily  the  result  of  formation  and  propagation 
of  four  types  of  rail  flaws 
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• Transverse  fissures  (Figure  3) 

• Horizontal  split-heads  (Figure  3) 

• Vertical  split-heads  (Figure  3) 

• Bolt  hole  cracks  (Figure  4). 

The  location  and  shape  of  these  flaws  and  recommended  actions  until  the 
flawed  rail  is  replaced  are  shown  in  Figure  5. 

The  transverse  fissure  is  a vertical  elliptical  defect  which  is 
often  totally  embedded  within  the  rail  head.  These  flaws  can  cause  complete 
rail  fracture  and,  therefore,  derailments.  They  are  particularly  difficult 
to  detect  because  thermally  induced  compressive  stresses  in  the  rail  compress 
the  flaw,  making  it  nearly  invisible  to  detection  techniques. 

Horizontal  split  head  defects  consist  initially  of  flaws  totally 
embedded  within  the  rail  head;  however,  these  flaws  "break  out"  of  the  head 
and  can  grow  to  3 to  4 feet  in  length.  At  this  later  stage  in  its  develop- 
ment, for  some  unknown  reason,  the  two  ends  of  the  crack  can  turn  up  and 
down,  respectively,  fracturing  the  rail.  This  also  results  in  derailment 
of  the  vehicle. 

Vertical  split  head  defects  are  vertical  separations  totally  with- 
in the  rail  head  which  causes  loss  of  running  surface. 

Bolt  hole  cracks  and  head  web  separations  arc  flaws  that  initiate 
at  bolt  holes  or  beneath  the  Lead  at  the  rail  end  at  a bolted  joint  con- 
nection. 

The  loading  cycles  that  these  flaws  experience  are  very  complex, 
in  particular,  in  the  head  region.  Both  vertical  and  lateral  loads  can  be 
simultaneously  applied  to  the  rail  and  these  loads  may  vary  greatly  from 
car  to  car.  The  stress  distribution,  Figure  6,  consists,  in  part,  of  a 
bending  stress  which  completely  reverses  during  the  loading  cycle,  and 
contact  stresses  which  are  highly  compressive  (-160,000  psi)  during  the 
loading  cycle.  The  transverse  shearing  stress  beneath  the  contact  region 
completely  reverses  during  the  loading  cycle.  The  problem  is  complicated 
by  plastic  flow  that  may  continue  over  a large  number  of  cycles.  Residual 
stresses  are  known  to  be  present  in  the  rail  head  and  have  been  measured  to 
continuously  increase  after  10^  cycles  (which  could  be  approximately  1/3 
year  of  service). 

It  is  thought  that  much  of  a rail's  life  is  spent  in  initiation  of 
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FIGURE  4. 
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of  these  flaws.  Once  cracks  have  developed,  it  is  believed  by  some  investi- 
gators that  the  crack  propagation  rate  is  significantly  affected  by  the  mean 
value  of  the  stress  field.  This  is  in  opposition  to  the  traditional  view 
that  fatigue  crack  growth  is  most  affected  by  the  value  of  K.  This  effect 
is  going  to  be  examined  in  ongoing  research. 

Naturally,  since  the  stress  states  involved  are  generally  three- 
dimensional,  TSC  is  interested  both  in  experimental  behavior  and  analytical 
solutions  to  mixed  mode  crack  propagatii  n.  It  was  also  remarked  that  the 
retardation  problem  is  also  of  interest. 

Oak  Ridge  National  Laboratory  (ERDA,  G.  Smith) 


Oak  Ridge  is  presently  investigating  reliability  of  flawed  pres- 
sure vessels  as  connected  with  light  water  reactor  safety. 

Reactor  vessels  are  in  general  200  inches  in  diameter  having  wall 
thickness  of  6 to  10  inches.  Operating  temperature  for  the  vessel  is 
approximately  500  F.  From  a fracture  point  of  view,  two  areas  of  interest 
are  at  nozzles  and  valves. 

Material  characteristics  for  the  vessels  are  difficult  to  obtain 
due  to  the  extremely  large  thickness.  In  order  to  get  values  even  a 
12th  specimen  is  inadequate. 

Oak  Ridge  has  performed  various  model  tests,  such  as  epoxy  models 
to  examine  flaws  at  inside  corners.  Numerous  model  pressure  vessels  have 
been  also  constructed.  In  tests  involving  pressure  to  failure  typically 
the  vessels  failed  at  3 times  the  design  pressure.  Flaws  that  were  modeled 
in  these  tests  were  typically  spherical  2 to  3 inch  flaws. 
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SESSION  VI 
BENCHMARK  PROBLEMS 


SUMMARY 


i 

I 


In  preparing  for  the  Workshop,  each  of  the  participants  were 
asked  to  suggest  possible  benchmark  problems.  A number  of  participants 
responded.  These  suggestions  were  roughly  in  three  categories:  (1)  sur- 

face cracks  in  a slab,  (2)  through  cracks  in  a place,  such  as  a CT  speci- 
men, (3)  corner  cracked  holes.  Task  forces  were  set  up, composed  of  interested 
participants, to  consider  each  of  these  problem  classes.  The  remaining  parti- 
cipants discussed  documentation  standards  for  reporting  benchmark  problem 
results.  The  resulting  documentation  standards  are  given  below.  Each  of 
the  task  groups  reported  on  their  deliberations  to  the  entire  group.  After 
some  discussion,  benchmark  problems  were  agreed  on  in  each  of  the  iree  cate- 
gories. These  consisted  of  problems  in  which  ranges  of  parameters  were  op- 
tional. However,  in  each  case,  one  specific  geometry  was  chosen  as  the 
mandatory  benchmark.  The  three  geometries  consist  of  the 

(1)  The  surface  crack  (circle  on  semi-ellipse)  in  a slab 

(2)  The  corner  cracked  hole 

(3)  The  compact  specimen. 

Specifications  for  these  problems  are  given  below. 


V^t^S’VS ■ 
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CRITERIA 

I.  Procedure 

1.  Describe  method  of  analysis 

2.  Describe  data  reduction  procedure 

3.  Give  estimate  of  accuracy 

II.  Validation  Steps 

1.  Tv?o  or  more  analyses  in  agreement  to  (?)% 

2.  Experimental  verification 

III.  Results  Documentation 


1. 

K 

Mandatory 

2. 

ui 

Mandatory 

3. 

Desirable 

4. 

Computer  effort  (in  cost) 

a.  Computer  cost  and  cost  algorithm 

b.  Manpower  required 

IV.  Full  Problem  Description 

1.  Boundary  conditions 

2.  Loads 

3.  Geometry 
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BENCHMARK  NUMBER  1 
THE  SURFACE  FLAW 


1.  Figure  1 shows  the  defined  germetric  nonmenclature 


la 


2.  Loading  shall  be  remote  tension 


a = a , z = K 
zz  o'1 


3.  Material  properties 


v = 0.3,  0.49 
E = 107  psi 


4.  Geometries 


A.  Semi  Circular  Crack 


~ = 0.5,  H/W  > 2,  \ < 0.2,  7 > 5 

A-C  u C 

B.  Part  Circular  Crack 

“ - 0.25,  ^ 0.70* 

C.  Semi  Elliptical  Flaw 


—■  = 0.25,  * * 0.25,  £ = 0.75 
2c  t t 


hi 


* This  dimension  was  proposed  as  0.75  at  the  Workshop.  However,  a number 
of  experimental  and  analytical  solutions  have  been  obtained  for  a/t  = 0.70. 
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5.  Presentation  of  Results 


In  accordance  with  documentation  standards,  plots  and  tables  should  be 
given,  wherever  possible  of  the  following  results: 


A.  in  10  increments  of  <Ji. 

o /ira 
o 


Note:  For  the  semi  elliptical  crack,  <j>  is  the  "eccentric  angle" 

measured  to  the  point  on  the  circle  of  radius  a with  the 
same  x-coordinate  as  the  point  on  the  ellipse. 


B.  C.O.D.  on  y = 0 and  x = 0. 


C.  Displacements  on  the  crack  surface  at  10  increments  of  <}>  on  con- 
centric ellipses  or  circles  with 


0.1  (0.1)  0.9 


a 0.95,  0.93 
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BENCHMARK  NUMBER  2 
THE  CORNER  CRACKED  HOLE 


1.  Figure  2 illustrates  the  nomenclature 

2.  Loading  shall  be  remote  tension 


a = o , z = H 
zz  o'1 


3.  Material  properties 
v = 0.3,  0.49 
E = 10^  psi 


4.  Geometries 


2R  + c) 

H - 2 tt 

a/c 

c/R 

a/t 

0.5 

0.2 

0.1 

1.0 

0.5 

0.5 

1.5 

1.0 

0.75 

2.0 

0.2 

0.2  (Mandatory) 

5.  Presentation  of  results 
K o 

A.  I fo:  10  increments  in 


(N.B.  K = S.I.F.  at  angle  <j>  for  an  imbedded  ellipse  in  an  infinite 

ivj 

00 

medium  for  o = a ) 
z o 

B.  Displacement  and  stresses  on  Lhe  crack  surface  at  10°  increments  in 
<j>  and  at  various  radii  position  (c.f.  Benchmark  Number  1) 


aty.n,. 


BENCHMARK  NUMBER  3 
THE  COMPACT  SPECIMEN 


Considerable  discussion  was  held  on  this  specimen  with  regard  to 
the  geometry  of  the  chevron  slot  machined  into  the  specimen.  In  conversa- 
tions with  C.  F.  Feddersen  of  Battelle,  who  is  Chairman  of  Subcommittee  I 
on  Fracture  Test  Methods  of  the  ASTM  E-24  Committee,  it  was  established  that 
the  ASTM  S399-74  standard  approach  to  preparing  a compact  specimen  is  to 
machine  a "convex"  chevron  notch  with  is  subsequently  fatigue  loaded  until 
the  fatigue  crack  has  propagated  into  a nearly  straight  through  crack  inter- 
secting the  surface  at  a distance  a from  the  applied  load.  Thus,  the  intended 
te;  t specimen  contains  a straight  through  crack  and  this  is  proposed  as  the 
standard  benchmark.  Possible  variations  on  this  problem  'lode  smoothly 
curved  cracks  to  account  for  variability  in  the  shape  of  the  crack  grown  by 
fatigue. 

Thus,  the  following  standard  benchmark  is  proposed: 

1.  Figure  3 shows  the  geometrical  nomenc] ature 

2.  Material:  v = 0.3,  0,49;  E = 10^  psi 

3.  Loading  P = 1000  lb/in,  through- the- thickness  point 

loading  at  x = -a  = =U.5W,  y = 0.275W 

4.  Geometries 

H = 0.6W,  t = W/2 

(other  dimensions  shown  in  Figure  3) 

5.  Presentation  of  results 

a.  Kj  as  a function  of  z,  0 z <_  t 

b.  Normal  displacements  and  stresses  for  y = 0;  z ~ 0,  t/4, 
t/2 ; -0.75W  < x < 0. 


FIGURE  3.  THE  COMPACT  SPECIMEN 


